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Abstract 

The full theoretical analysis of the kinetics of multicomponent nu- 
cleation is presented. The relief of the free energy with surface excesses 
was analyzed, the valleys and ridges were described, their mutual in- 
teraction was studied. The new possibility to change the valley of 
nucleation is shown. The possibility to have one common valley in- 
stead of several neighbor ones which leads to the radical change in the 
height of the effective activation barrier and to the new value for the 
nucleation rate. 



Introduction 

Historically the problem of determination of the stationary rate of nucleation 
was primary investigated in one- dimensional approximation [1], [2]. The 
unique variable characterizing an embryo of a new phase was a number of 
molecules inside the embryo. Meanwhile, it is evident that the embryo has at 
least several characteristics, which have to be taken into account to give the 
adequate description of the nucleation process. That's why it is necessary 
to study the description of the nucleation process on the base of several 
characteristics of the embryo. 

One can not pretend to take into account all characteristics of an embryo 
and to give the nucleation description on the base of all embryo characteris- 
tics. The kinetic aspects of the embryo formation are also far from clear inter- 
pretation. Even the mathematical structure of the theory of multi-dimension 
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nucleation is far from complete understanding. So, it is worth to start with 
the simple cases of multidimension description. 

The simplest and the most evident example of mult i- dimension descrip- 
tion is a multicomponent nucleation. It means that the nucleation in mixture 
of vapors is studied. Here kinetic coefficients are determined extremely clear, 
the free energy of the embryo is also rather well known in general features. 

The history of investigations of the binary nucleation is very rich. The 
number of publications concerning the binary nucleation is now greater than 
devoted to other domains of the nucleation theory. But already in the case 
of binary nucleation there appear many problems to solve. So, it is worth 
paying attention namely to multicomponent nucleation. 

Until nowadays there is still no universal self consistent analytic approach 
which makes use of all previous theories or directly shows their errors. This 
task will be the goal of the present paper. 

At first one has at least to mention approaches, which pretend to give 
original recipes for the stationary nucleation rate. Certainly, the classical 
expression for the free energy given by the standard thermodynamics has 
to be the starting point of a theory. In our analysis we ignore approaches 
suggesting some artificial correction terms or some reconsiderations without 
a solid thermodynamic base. 

The microscopic corrections to the free energy given by classical thermo- 
dynamiocs [5], [1] are not the subject of our investigation, we consider only a 
task to describe nucleation at the relatively low supersaturations. Even this 
question is out of a true solution. We do not consider a normalising factor 
in an equilibrium distribution which evidently appear in the expression for 
the nucleation rate. This will be a subject of a separate investigation. 

The first essential contribution to establish the binary nucleation rate 
was made by H.Reiss [5] who determined the rate of nucleation on the base 
of a steepest descent line in a near critical region. Solution of a kinetic 
equation presented by Reiss was corrected by Stauffer [6]. In the last paper 
the correct formula for the rate of nucleation in the square approximation of 
the free energy in the neighborhood of the critical embryo was given. Earlier 
the general ideas for the problem of overcoming the activation barrier in the 
multicomponent case were formulated by Langer [7j but one can not state 
that the publication [6] is a direct consequence of [7]. One has to stress that 
the constant direction of a fiow in a neighborhood of the critical point was 
simply postulated in [6]. This constancy can be proven only with the help of 
the boundary conditions which was done in [S] where direct solution of the 
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kinetic equation was presented. 

As it became clear after the solution of Trinkaus [H] the problem to deter- 
mine the nucleation rate requires to decide whether the transition over the 
barrier really occurs at the critical point ( the saddle point - here and later 
the critical point means the coordinates of the critical embryo). When there 
is a strong hierarchy between kinetic coefficients of absorption of different 
components one can see that flow of embryos can pass aside the critical re- 
gion (the region near the critical point), but over the ridge far from critical 
point. 

Solution of Stauffer implies the square approximation of the embryos 
free energy near the critical point. Solution of Trinkaus implies the linear 
approximation of the height of special activation barrier. But as it will be seen 
later there is no contradiction between approximations - both are suitable in 
corresponding situations. 

Further analytical progress is associated with the appearance of many 
variations. Among them one can outline the refined Stauffer's solution pre- 
sented by Berezhkovski and Zitserman [TU] and conception of the genuine 
saddle point proposed by Li et al. [TT]. One has to stress that these contri- 
butions did not radically change the already known formulas for nucleation 
rate but slightly corrected some known results. In this context it is also 
necessary to mention the publications of Shi and Seinfeld [12] and Wu [13] . 

Here we do not analyze the theories connected with the reconsidered 
free energy of the embryos formation taking their history from the famous 
publication of Lothe and Pound |^14j| and modifications of this approach to 
the case of multicomponent nucleation. Any new expression for the free 
energy will cause the new value of the nucleation rate but the mathematical 
structure of the derivation of the nucleation rate remains the same. 

In the middle of 1990-ies the serious set of attempts to analyze the binary 
nucleation problems was presented in [15], [T6] . 

In the last years one can outline the publications which analyze the same 
problems which have been already mentioned. The problem of boundary 
conditions was revised in publication of Wilemski and Fisenko [T7]. The 
authors put the natural boundary conditions directly at the boundaries of 
a whole pre-critical region where these conditions are evident. But then it 
is necessary to solve the kinetic equation in the whole pre-critical and the 
near-critical region which was done in [13| only numerically. 

The set of papers by Li, Nishioka, Maksimov [11] is devoted to give the 
definition of the generalized saddle point which can be used both in the 
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situation where the flow goes over the standard saddle point and in the 
situation of hierarchy where the solution of Trinkaus |9] takes place. This idea 
is certainly attractive but as it will be shown in this publication sometimes 
the nucleation occurs in a more complex way and can not be described in 
terms of the genuine saddle point even approximately. Moreover, the point 
of the Trinkaus' solution depends not only on the free energy but also on the 
derivatives along special directions. 

The problem of transition of the binary case to the unary one was studied 
in [18] where the full analysis of this problem was given. Here we are not 
interested in this transition because the embryo with one molecule of a rare 
component can not be considered on the base of a standard thermodynam- 
ics in an approximation of a homogeneous liquid which is adopted in this 
publication. 

Here we do not analyze numerous publications which combine the stan- 
dard known approach with some artificial additions. Such combinations are 
rather typical for publications of Djikaev with coauthors (see, for example, 
|19j). In [in] the values of kinetic coefficients from the first passage time 
analysis are formally injected in the standard solution presented in [20] and 
the final formulas are presented. One has to mention that the first passage 
time analysis is based on some unknown characteristics (for example, the 
height of activation barrier for a molecule to penetrate inside the embryo) 
which can not lead to concrete results. 

All mentioned publications in the binary nucleation do not make any 
profit from the topology of the relief of the free energy of the embryo. This 
task was solved in ^20j where the structure of relief of the free energy in 
the capillary approximation was studied. It was shown that the relief of the 
free energy can be characterized as the straight channels, ridges and saddle 
points. In this publication the results of [20j will be widely used. 

The formulation of the capillary approximation faces the difficulty known 
as the Renninger-Wilemski's paradox [21], [22]. Because of publications by 
Oxtoby and Kashiev [23] the thermodynamic background of the surface ex- 
cesses is completely studied. To overcome this difficulty one has to write the 
Gibbs' absorption equation and to introduce surface excesses of components 
of the surface of tension. This leads to the difference of concentration in the 
surface layer and in the bulk of the embryo which was noted in [20]. But 
there further conclusions for kinetics of the process have not being made. 

The structure of the free energy relief with surface excesses was investi- 
gated in pS] but only in thermodynamic aspects of the problem. The kinetic 
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features have not been considered in [23] . 

In [8] the kinetic equation was solved in the neighborhood of the critical 
embryo. The progress achieved in [8] was the appropriate formulation and 
account of boundary conditions. Certainly earlier the boundary conditions 
were mentioned in [9] but they were put in the infinitely far points where the 
structure of a free energy can not be seen in all details. Namely the necessity 
to conserve the boundary conditions at the low boundary of a near-critical 
region determines the conservation of the square form of the free energy in 
transformation presented in [8j. 

Having summarized the development of theoretical investigations in the 
binary nucleation one can state that despite the essential progress in this 
field there are still many problems to consider. 

It is rather natural to construct the global picture of the nucleation in- 
cluding the case of the hierarchy between kinetic coefficients, surface excesses, 
etc. The unification of the free energy topological features with the already 
mentioned approaches is the main goal of this publication. This approach 
leads to many rather essential features of nucleation presented below. More- 
over, some striking features changing the rate of nucleation in the order of 
magnitude will appear. 

The structure of this paper is the following 

• In the first part the main ideas of the capillary approximation are 
formulated and the free energy is constructed. Here the surface excesses 
are taken into account. The variables providing the simple form of 
the free energy are shown and their connection with the numbers of 
molecules in the embryo is established. 

• The second part is devoted to the description of the near-critical region. 
It is shown that this region has the form similar to the case of the 
absence of the surface excesses. Here the hierarchy of evolution will be 
shown. 

• The third part is devoted to the analysis of the Reiss' solution and 
the Stauffer's one. The plausible way to see the Reiss' formula will be 
shown. The moderate value of the difference between the Reiss' and 
the Stauffer's solution is justified. This is important for the possible 
ignorance of the slow or rapid variables of correcting order. 

• The forth section analyzes the jump of the embryos from one channel 
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to another one. The situation of the near-equihbrium faUing transition 
is considered here. The solution is found also in the discrete model. 

• The fifth section considers the conception of the common valley. The 
equilibrium common valley transition is analyzed. It will be shown the 
new height of activation barrier. This value seriously differs from all 
known results. 

• The sixth section analyzes the general picture picture of the nuclcation 
rate formation. The case of the equilibrium saturation of the destina- 
tion valley is studied. 

• All results are summarized in the conclusion. 

1 Thermodynamic basis 
1.1 Capillary approximation 

The main object involved in determination of the nucleation rate is the free 
energy of the isolated embryo. To give the description of the embryo one 
has to fix the variables of the state of the embryo. Assuming the thermal 
equilibrium of the embryo one can describe the embryo only by the numbers 
of molecules Ui inside the embryo. These variables are extracted by following 
properties 

• In elementary acts of evolution Ui are changed separately. The step of 
change is one unit. 

• Although the free energy even in capillary approximation is not diago- 
nal the form of expression for the free energy is relatively simple. 

In the capillary approximation the energy F ordinary taken in the thermal 
units is the sum of the bulk part B and the surface part fl. 

F^-B + Q 

The ordinary expressions for B and Q are following 
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Here the sum is taken over all components of the embryo, /li are the differ- 
ences of the chemical potentials counted from the equilibrium values (with a 
negative sign), S is the square of the surface of tension, 7 is the renormalized 
surface tension. 

The difference between the precise value of the free energy and the value in 
the capillary approximation referred as "correction terms" (c.t.) is supposed 
to be relatively small in comparison with B + This situation takes place 
when the number of molecules 

J^tot = X] ^» 

i 

inside the embryo is very (strictly speaking infinitely) big 

l^tot > 1 (1) 

The inverse number of molecules (i.e. Vtot) will be the small parameter of 
the theory. So, 

F = -B + n + c.t. 
where c.t. indicates correction terms with a property 

\c.t.\ < |F| 

Ordinary the decomposition of correction terms on inverse radius of the 
embryo converges and F has the form 

00 

F ^ -B + n+ CkT'^ + Co In r 

k=-l 

Here Ck are the coefficients. One can also consider the last decomposition as 
an asymptotic decomposition. We shall accept the validity of this decompo- 
sition. 

Ordinary this decomposition is taken with a finite number of terms 

lH-2 

F = -B + n+ Ckr-^ + colnr (2) 
fe=-i 

Prom the last decomposition it follows 

\d c.t. I <C \dF\ 
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The last inequality is important in the justification of the linearization of the 
free energy. 

While speaking about the capillary approximation one has to imply a 
whole set of assumptions beside the pure thermodynamic consideration. There 
are several ordinary used approximations included into the capillary approx- 
imation. These approximations are the following 

• The surface tension is attributed to the dividing surface calculated on 
the base of the volume separation, i.e. 

i 

where Vi are the volumes in a liquid phase. The formal factor An/ (47r/3)^/^ 
is ordinary included into the effective surface tension. 

• Values Vi, 7 are taken from the case of a bulk liquid. 

• To give expressions for Hi one has to use some model. The most widely 
used model is the model of a liquid solution. The validity of this model 
requires 

> 1 

for every component. Certainly one can use other models and ignore 
these limitations. When uj = 1 for some component, one can consider 
this component as a heterogeneous center. That's why the extension 
of the approximation of the regular solution up to z/j = 1 in JTE\ causes 
questions. 

In this paper we shall use the model of solution. 

To give a formula for chemical potential one can define a supersatura- 
tion as 



where rij is the molecular number density in the existing vapor and rijoo 
is the molecular number density of the vapor saturated over the pure 
bulk liquid of component i with a flat surface. Then 

/i, = ln(C)-lne.-ln/,({0) 

Here it is supposed that the vapor is an ideal gas which gives the value 
for the first term in the r.h.s. as ln((^j). Ordinary it is assumed that the 
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concentrations form a set {^} of concentrations and the coefficients 
of activity fi can depend on the whole set of activities. 

For approximation of ideal solution all coefficients 

/i = l 

To know fi one has to construct some model of solution or to use some 
experimental data. 

A special question concerns the definition of concentration. Ordinary 
the concentration is determined as 

^^ = ^ (3) 

Sometimes this definition is also included into the ordinary auxiliary 
approximations of the capillary approach. 

This question is directly linked with the Wilemski-Renninger's paradox 



All assumptions made above are necessary for formula for the free energy 
in capillary approximation. 

The presented formula for the free energy is rather transparent, but it 
faces the difficulty known as the Wilemski-Renninger's paradox. The diffi- 
culty is the following: 

• It is known that in the critical embryo the Kelvin's relation 

Mi 

— = invariant 

has to be observed. This follows from the general thermodynamics and 
from the sense of chemical potentials. 

The last relation gives an equation for the concentration in the critical 
embryo. 

• One can come to the same equation on concentration also by direct 
differentiation of expression for the free energy. For simplicity assume 
that Vi do not depend on concentration. This gives 

dF 

dii 
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Here it is supposed that Vi do not depend on concentration. 

One has to recall that the coefficients of activity fi satisfy the Gibbs- 
Duhem's equations 

^ iidni = 

i 

which put a restriction on the coefficients of activity. 

i 

Moreover the Gibbs-Duhem equation can be written as 

Then the differentiation becomes very simple and leads to 

Then one can come to the widely known Kelvin's equation only if the 
derivative d'~f/dC,i is zero. So, the formal recipe is to forbid the differen- 
tiation of the surface tension on concentration. Since the last equation 
comes from the foundations of thermodynamics it means that some- 
thing is irrelevant in the previous formula for the free energy. 

In the case when Vj depend on concentration we have 

dF sr^dnj 2 x-i/sr sr^dvj , ^d'y 

8^ = = - E + 57(1: '"l". + E + 

But the Gibbs-Duhem equation has to be here the following one 

_ y + 2 V ^ ^ )-i/3 V = 

which leads to the same conclusions. 
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As the result there appeared a formal recipe not to differentiate the 
surface tension. At first it was the artificial recipe but later the jus- 
tification of this recipe was given on the base of the Gibbs dividing 
surfaces formalism. 

To resolve this difficulty one has to add to the free energy some new 
contributions connected with the surface excesses. It will be done later. 

Now we return to consideration of the properties of F. 

The leading idea here is the extraction of the mentioned small parame- 
ters Recall that conditions are rather far from the second order phase 
transition. If we accept that the surface layer has a finite thickness d, then 
in the limit r — > oo (where r is the radius of the embryo) one can see that 
correction terms {c.t.) are really relatively small in the following sense 

\B\ > \c.t.\ |Q| > \c.t.\ 

I dc.t. I ^ I I I ^c.t. I ^ I I 
dvi dvi dvi di^i 



These inequalities are valid for absolute values. 

These inequalities is a new result and they will be widely used below. 
Their validity can be proven analytically. 

1.2 The form of the free energy 

To see the structure of the free energy one can introduce the extensive variable 

i 

Certainly ^ is a volume of the embryo. Then 
with the generalized chemical potential 
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The generalized chemical potential allows an interpretation 

< V > 

as the ratio of the mean chemical potential excess and the mean volume per 
one molecule in the embryo. 

One can take also as an extensive variable the total number of molecules 
inside the embryo 

i 

Then the free energy can be written in a following way 

i 

where 

i 

and the renormalized surface tension 

i 

appears. 

One has also to mention the possibility to take as external variable the 
surface energy in the power 3/2, i.e. 



used in Then the free energy has the form 
where the generalized chemical potential is 



The third possibility used in [20] is the most preferable because here the 
free energy has the most simple form and in the "surface" term the factor 
depending on concentration is absent. 
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The question to discuss is what we shall take as an extensive variable 
- the variable proportional to the volume or the value proportional to the 
number of molecules? Thermodynamics does not give an answer because 
asymptotically these values are proportional. 

But the problem to take into account the Renninger-Wilemski's paradox 
remains here our of attention. To overcome this difficulty one has to include 
into description the surface excesses of components. To take these excesses 
one has to choose the surface accurately. The most preferable choice is to 
choose as the surface the surface of tension because the surface tension can 
be attributed to this surface without corrections. At this surface all compo- 
nents have the surface excesses ipi but the surface tension can be attributed 
namely to this surface. In the first (rough) approximation these values are 
proportional to the square S of the surface of tension 

ipi = QiS 

Parameters Qi are supposed to be independent on S and have to be given 
by the theory of a liquid state. 

The square of the surface of tension can be approximately calculated as 

i 

Certainly, there exists a difference between a surface of tension and the sur- 
face covering the volume of the embryo. But since the sense has only 5*7 one 
can attribute this difference to the value of 7. 

Here we omit the constant factor having included it into the surface ten- 
sion 7. Hence, 

i 

The last relation is not a formula for ipi but an equation. It can be solved 
by iterations. These iterations are based on a small parameter ipi/vi. The 
smallness of these parameters at z/j — > cxo is evident. The first approximation 

i 

is already suitable as a leading term under the conditions ([1]). The second 
iteration 

i j 
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will refine the solution. The complexity of dependence of ijji on is the 
certain difficulty. 

The value of concentration has now to be redefined as 

^ _ ^i- QiS 

si ~ 



As an extensive variable it is natural to choose the straight analog of q, 
namely 

i 

But this choice does not lead to the "true" form of the free but to 
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with 

Aj = i/i - V'i 

Here the dependence 7 on {^} is the source of difficulties. Certainly, 

— ^ 

One can introduce another set of variables. Now instead of ^ one has to 
choose the extensive variable 



In these variables the free energy F has the form 

F = -KbgiO + (4) 
with the generalized chemical potential 

K 

or 

i ^ 
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One has to show that bg does not depend on k. To fulfill this derivation one 
can come to 

"9 ~ Z^^»/^« C3727I 



or 



It can be also presented as 

The last relation evidently shows that bg is really a function of ^. The 
dependence on k is absent. 

One can use expression ([5]) to clarify the Renninger-Wilemski's paradox. 
According to the Gibbs' absorption relation 

d'j = dJ2 Qil^i 
j 

the derivative of the surface tension on concentration is cancelled by the 
corresponding derivatives of qi on ^. So, if we write bg without surface 
excesses as ^ 

we have to forbid the differentiation of 7 on concentration. Now the Renninger- 
Wilemski's paradox is explained. It is necessary to stress that the reason is 
not the formal Gibbs' absorption equation, but the difference of concentra- 
tions in the bulk solution from the integral values. 

Although the the new variables ensure the simple form of the free energy 
their connection with "initial" variables Vi is rather complex. One has to see 
how on the base it is possible to reconstruct z/j. The procedure is the 
following: 

• On the base of we know /i^, then we get 7 — I]i QiHi- 

• This gives a value of 

S = «:2/3/(^ _ ^ ^^^^) 
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• On the base of S having presented S as 

S = J2 ^i^i = J2 J2 

i i j 

we get J2i Ai- 

• Since Aj = J2j we get all Aj. 

• Then 

Ui = \i + Qii^S 

and we know all z/j. 

The inverse transformation can not be made by explicit formulas, the 
problem to fincQ g on the base of u has been considered above. When g is 
found then A is known. This gives ^ and k. 

The main new facts found here are the following: 

• The variables giving the simple expression for the free energy with 
surface excesses are found. 

• The recipe to get the initial variables on the base of the new ones is 
given 



1.3 The structure of the free energy reUef 

The functional form (jlj) has some consequences analogous to those considered 
|20j . But now this form takes into account the surface excesses of an 
embryo. Here the form (jl]) ensures the following properties of the free energy 
of an embryo 

• One can see the channels of nucleation defined by equations 



< 



Along these channels the equilibrium density of distribution has a max- 
imum (but the real distribution coincides with the equilibrium one only 
in the part of the pre-critical region) 



^When the index is absent it means that the whole set is considered. 
^In [20] the free energy without surface excesses was considered. 



16 



• Because of the Gibbs-Duhem's equation the variables in differentiat- 
ing of bg are separated. This leads to the approximately zero value of 
the cross derivatives d%g/ d^id^j. 

• One can see the separation lines of nucleation defined by equations 
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Along the separation lines the equilibrium density of distribution has 
a minimum. 

In one channel there is only one saddle point. Certainly, this takes place 
only in the capillary approximation. This saddle point has a coordinate 
Kc determined from the following equation 

o 



'^bg{Q 

Here is the coordinate of the channel. 

• The amplitude value of the free energy Fc in the channel is given by 
the formula 

F = -K^/^ 
" 3 ^ 

Here one can see the Gibbs' equation and now it is clear that namely 
K^/^ is the true surface energy, but not as it seems from the first 
point of view. One has to attribute to the surface energy all energy like 
contributions with the space dimension 2 (or 2/3 in relative units). 

• All channels are independent - the embryos starting from the origin of 
coordinates will use only one separate channel to go to the supercritical 
region where they begin to grow irreversibly. The nucleation flow will 
mainly go through the channel with minimal Kc or maximal bg. This 
remark concerns the case where there is no strong hierarchy between 
kinetic coefficients of absorption. 

One can see that the picture of nucleation is rather simple, but this sim- 
plicity was observed for the free energy with the surface excesses for the first 
time here. This is the new result of this section. 
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1.4 The form of the near-critical region 

As it has been mentioned at the beginning the set of natural variables is Ui. 
The elementary kinetic act of absorption leads to the change 

So, it is necessary to establish connection between k, ^ and 1/ at least approx- 
imately. 

Denote by the subscript o the values when all surface excesses are zero. 
Then the theory is very simple and one can get the connection between Kg, Co 
and 1/0 in a very transparent manner. Prom uq to Kq, one can get by 

1^^^ = 7 Vii'io 

i 

j 

Inverse transformation is given by the chain formulated above. So, it is quite 
easy to write the kinetic equation for the case of the absence of excesses. 

The above consideration shows the role of the case with zero excesses. 
Hence, this case will be the base to construct the description in the general 
case. 

Return now to the general case. 

One can define the near-critical region as the region where 

|i^-i^c| < 1 

This is quite analogous to the one component case. But here we consider the 
near-critical region associated with the given channel. Then it is necessary 
that this region has to be closer to this point than the separation lines. 

One can define the positive size Ak of the near-critical region along the 
channel of nucleation as 

Here is the coordinate of some channel. Certainly, we get two values AiK 
and A2K. corresponding to the positive and to the negative shift. In the 
square approximation of the free energy 

AiK = A2K = 3^2/3 (6) 
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When z/j ^ 1 for all i the square approximation is rather accurate. 

Analogously one can define the characteristic sizes A^j according to rela- 
tion 

Certainly, we get two values Ai^j and A2^j corresponding to the positive and 
to the negative shift. In the square approximation 

Aie. = A.e. = ||||-^/^«-^/^ (7) 

When z/j ^ 1 for all i and there is no singular behavior of generalized chemical 
potential then the square approximation is valid. 

Certainly, it is necessary that the channels have to be separated, i.e. the 
height of the separation line has to be several thermal units higher than the 
height of the channels. This has to take place at k, near the critical value. 

Wc define the reduced near-critical region as the region where |k — Kc| < 
An, — ^ic\ < A^i. This definition differs from the ordinary definition of 
the near critical region as extracted by condition \F — Fc\ < 1. 

In the multi-dimensional case there exists long tails near fines F — Fc. To 
illustrate it one can use the square approximation, then the curves F = Fc-|-1 
and F = Fc — 1 are hyperbolic ones with common asymptotics which are 
straight lines. 

Wc shall define the tails as the regions corresponding to \F — Fc\ < 1 and 
16 - 6c| > A^i, \k- Kc\ > Ak,. 

Actually, the following statements can be proven analytically: 

• One can show that the tails do not play any essential role in formation 
on the total nucleation flow. 

• Then it is possible to reduce the near-critical region up to the following 
domain 

\k — Kc\ < Ak 

Here and later we shall imagine the reduced near-critical region speaking 
about the near-critical region. 

Now one can see that the relative sizes of the near-critical region are small 
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Ordinary this smallness is implied when the kinetic coefficient of absorp- 
tion is supposed to be a constant value. Here this smallness will help to prove 
the following main result of this section: 

• In the near-critical region the function F — Fc as a function 

of variables k — Kc, ~ for every i has practically the same 
behavior as the function Fo — Fco as a function of variables 
iio — ^ico- At least the relative difference is small: 

The explanation and the idea of the proof is rather simple. Really, the 
correction terms to which the excesses belong begin to be essential only 
when the surface term cancels the bulk term. But as it clear from the 
sequential differentiation this can take place only in the first derivative. 
Starting from the second derivative the contribution from the bulk term 
is zero and this compensation can not take place. This effect is taken 
into account by the shift of instead of i^jco- So, here the influence of 
correction termms is negligible. 

The last result allows to write the kinetic equation in variables taking 
into account the surface excesses by a simple shift. This takes place only in 
the near-critical region. This result is new. 

1.5 The place of the Renninger-Wilemski's effect 

The " paradox" of Wilemski and Renninger occupies so important place in the 
multicomponcnt nucleation that from the first point of view it seems that this 
is the real effect taking place in the leading term of capillary approximation. 
Below it is shown that this effect has an order of correction. To see this effect 
one can redefine k as S^/'^ and forget about excesses. 
Really, from equation 

q^\k= fixed Q^l^ 
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it is seen that the first term with the derivative || has a correction 

order k^/^ 



in comparison with the second term -^k'j^^'^ having the order k 

We extract this result which is exphcitly outhned for the first time here 
because of its importance for the reconstruction of the logical self-consistency 
of thermodynamics. Only the correcting order of the term with the derivative 
of the surface tension allows to ignore it in the main order and to return the 
leading role of the ordinary capillary approximation. 

Since the formal recipe to resolve the Rcnninger-Wilemski's paradox is to 
forbid the differentiation of 7 on concentration then the equation on concen- 
tration will be different. It would cause the impression that there is a shift 
in a leading term. The correct answer is that this result causes the shift in 
Fc which has a correction order as it follows from the last equation. 

The necessity to develop the theory with surface excesses is evident be- 
cause the surface excesses will essentially shift the position of the near-critical 
region. The shift is many times greater than the size of the near-critical 
region. The shift has the order k (because there is another equation on 
concentration - the derivative of 7 on ^ is cancelled) while the size of the 
near-critical region has the order k^^"^. 

One can treat the surface tension as a coefficient in the first correc- 
tion term proportional to the surface of the embryo. The coefficients at 
K^/^, luK, etc. depend on intensive variables (concentrations is one 

example). Their derivatives will be cancelled by derivatives of correspond- 
ing excesses. The structure will resemble the Rcnninger-Wilemski's paradox. 
But here the dimension of "surface" will be K^^^,\nK, k^^^^, etc. This effect 
will be called as "generalized cancellation of derivatives on intensive vari- 
ables". 

One to note that the same procedure can be effectively apphed for all 

other correction terms. Rigorously speaking to determine the form of the 
near-critical region one has to take the expression for F with correction terms 
up to the order which causes the shift of position of the near-critical region. 
Now it is clear that the effect of all correction terms will be quite similar to 
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the already described one. 



2 Channels and separation lines 

2.1 Similarity of the near-critical relief 

Although the Renninger-Wilemski's effect has a correction order it is worth 
taking it into account. The main reason is the following: 

• The relative sizes of the near-critical region is very small. Really, from 
([2D it follows that 

— ~ n-'/' « 1 
and the relative size in K-scale is small. From ([7]) it follows that 

and the relative size in the C,i scale is small also. Then it is clear that 
the relative size in j/j scale will be 

/ -1/3 

Namely these estimates allow to put in the near-critical region the 
kinetic coefficient of absorption of the molecule of i-th component 
to the constant value Wj^ corresponding to the critical embryo 
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So, the relatively small error in the determination of the coordinates I'i can 
remove embryo out of the near-critical region which makes the consideration 
of kinetic equation without surface excesses in the near-critical region useless. 

Beside this one has to take into account that the elementary transitions 
are written in the u-scale 

Vi — * Z/j ± 1 

and the free energy is written explicitly (with the surface excesses account) 
in variables k, ^j. So, it is necessary to have the a very precise transformation 
between z/j and This forms the problem. 

Although the transformation from K,S,i to z/j exists it is very complex. 
The inverse transformation has not been found explicitly. So, it is necessary 
to establish the approximate connection. The following statement establishes 
this connection 
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• The function F — as a function of z/j — z/j^ approximately coincides 
in the near-critical region with the behavior of Fq — Fq^ as a function 
of Ui - Uico: 

Fo{uicQ + Hi) - FcQ ^ F{uic + Vi) - Fc 

This property can be called as the approximate similarity of the free 
energy relief. 

Here this fact is established for {z/j} variables while earlier the same con- 
clusion was made for n, {^} variables. 

The idea of the proof of this property is based on the simple remark 
that the terms produced by the surface excesses can be important only when 
the terms produced by the bulk and surface contributions are cancelled. In 
the near-critical region this occurs only in the first derivative over i/^ at the 
critical embryo. Cancellation in high derivatives is impossibl^. Here we use 
the form F = J2if^i^i ~ iS and differentiate it over z/j. This ensures the 
similarity of relief. 

Now one can propose the following sequence of actions 

• at first one has to solve equations for the characteristics of the critical 
embryo 

• then one has to solve the kinetic equation without excesses but in 
shifted coordinates. 

Certainly, the similarity of relief takes place both in i/j coordinates and 
in K, coordinates. 

Analogously one can one can prove the small relative role of microscopic 
corrections in the value of dFc/ dQ which is used in construction of the global 
evolution of the phase transition. Here the formula ([2]) has to be used and it 
has to be taken into account that the coefficients in this formula are constant. 

2.2 The form of pre-critical region 

One can see the following important property: 

• The critical embryo can not have = 

■'Since the high derivatives of the bulk contribution are zero. 
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It can be seen from the explicit form of /Zj 



fii = lnCi + ln^i + ln/i(0 
Recall that here Q is the supersaturation of i-th component defined as 

Q = — 

Hi is the molecular number density of vapor of i-th component, na is the 
molecular number of the pure saturated vapor of i-th component. The sec- 
ond term is caused by the standard entropy of mixing, the third term char- 
acterizes the deviation of mixture from the ideal solution, here fi is the 
phenomenological coefficient of activity. 

Then one can see that at ^ 1 the situation of dilute solution takes 
place. Then the Henry's law states that the situation is close to the ideal 
solution, then fi = l and there are no correction terms. Then one can see 
that 

db„ 

and the condensation into the pure component is forbidden. Analogously 
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Earlier the analogous estimates were formulated in ^20j for C,io- Then from 
IQ and ([7]) it follows that the widths Aui along Ui satisfy 

Aui > 1 

for all i. These estimates ensure the possibility of continuous description of 
evolution in the kinetic equation. 

In the absence of the strong hierarchy between coefficients of absorption 
one can define the pre-critical region by two conditions 

• by inequality 

F < F — 1 

^ ^ ^ cm ^ 

where Fcm is the minimal activation barrier among different channels. 

• by requirement that this region has to be continuous and the origin 
belongs to this region. 
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One can prove that in this region the quasi stationary equihbrium state 
takes place. Here the absence of the hierarchy of kinetic coefficients plays 
the principal role. 

Now one can investigate the form of the pre-critical region in i/j variables. 
It looks like a star and the needles are going along the bottoms of channels. 
Certainly due to restrictions i^i > there is only one quarter of a star. In 
K, {^} variables it looks like a brush. 

If in every channel we put the value Fc corresponding to this channel, 
the shortest needle is the main one. The shortest needle (in k, ^ plane) 
corresponds to the lowest barrier and, hence, it is the main needle through 
which the nucleation takes place. 

If the level F^ is chosen as F^m and it is one and the same for all channels 
then the main needle is the longest one 

To see the relaxation to the equilibrium distribution we need to determine 
the minimal diameter of this star. It is given by the following relation 

l^minbg min ~l~ ^min — Pcm 1 

Here hg min is the minimal value of hg. So, if \hg min\ does not go to infinity, 
one can easily see the finite value of Kmin and the connection of channels. 

The last consideration solves the problem of connection of channels of 
nucleation. The problem was that the behaviour of channels near the origin 
where the surface excesses can play the leading role was unclear. So, one 
could not say whether the channels are connected or no. Now the concrete 
position of channels near origin is not important. 

The only condition is the restriction on - this function can not go to 
— oo at some concentrations. 



2.3 Characteristic sizes of near-critical region. 

Consider the variables parallel to ^i,K and having the scale of 1/^. These 
variables are 



7 / ^ 

^par 



^i perp — ^par^i 

Here it is supposed that all Vi have the same order of values. The total 
number of components n is not supposed to be a big parameter. 
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Then the halfwidths along Upar and Ui p^rp satisfy the following estimates 

^^par ~ ^ ^l^i perp 

AUpar ~ K^/^ ~ l^tlt 
A 1/2 1/2 

The time of establishing of the stationary state along Upar, perp is given 

by 

W+ A-l 



tr iperp ~ 



w+ , 

" '^i perp 



Here all kinetic coefficients of absorption are supposed to have one and the 
same order of value which is marked by . 

Then we come to the following strong inequality 

Upa^ _ ^1/3 ^ ^ 
tr iperp 

This equation states the hierarchy in the near critical region. Earlier this 
hierarchy was established in [25] for the situation without surface excesses. 
Here it is done for the presence of the surface excesses. 

To see this property the main effort was spent to show the similarity of 
forms of the free energy relief. Then one can come to the hierarchy rather 
automatically. 

The mean characteristic time tu to overcome the near-critical region for 
the embryo at the bottom the channel at the boundary of the near-critical 
and pre-critical regions has the order of tr par 

^r par 

Then we come to the following strong inequality 

ill 



tr iperp 



k'/' » 1 



It means that along z/j perp there is a quasi equilibrium. 
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2.4 Advantages of hierarchy 

On the base of hierarchical inequahties one can see that along Vi p^rp or 
there is quasi equilibrium. Then the distribution function n{{i'i}) which can 
be transformed into n{vpar-i {I'l par}) can be presented as 

where Npar plays the role of the amplitude of the known equilibrium distri- 
bution and rieq is given by 

^eg{{^i perp}) ^ exp( -F^l^pari {^i perp}')) 

or more convenient 

^eq{{^i perp}) 6Xp(i^(z/par ) {^i perp}) F (^Upar t {yi perp b})) 

where b marks the coordinate of the bottom of the channel. 

Then there remains only the task to determine the amplitude Npar- This 
is a simple onc-dimensional problem of nucleation. One can easily solve it. 

Reduction the problem of nucleation to the one dimensional case allows 
to solve more complex situations. At first one can see that when the char- 
acteristic width of equilibrium distribution seriously changes it leads to the 
change of the effective free energy in Npar- Really the effective free energy 
looks like 

F,ff -In AeqU 

where 

^eq^ — ( ^ ^ '^eqiyperp)) 

In the majority of cases the summation in the last formula can be replaced 
by integration 

^eq (^perp) dfperp) 

-oo 

Here the region of integration is formally put to an infinite one, actually one 
has to integrate over the region near the bottom of the channel where Ugq is 
essential. 

The further simplification is the following: one can take the last integral 
in square approximation for the equilibrium distribution: 

/J T\ / p TT ( ^ i^par 1 {^i perp}) I / \2\ 

neq[\T^i perpf) ~ CXpi^ r},) J_|expi^ r,o 2 perp=i'i perp b\^j perp perp b) ) 

j j perp 
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This allows to take integrals explicitly. 

Then the effective free energy is given by 



eff 



F-^ln 



IdlP^ perp — perp b 

J perp 



Later one has to solve one-dimensional nucleation problem with the effective 
free energy instead of the initial free energy. As it will be shown later by 
demonstration of the plausible derivation of Reiss' formula one has to be very 
attentive at this step. 

Ordinary in the near-critical region the value Aggi^ is constant and there 
is no peculiarities in behavior of i^e//- Certainly, in the square approximation 
the is an explicit solution of Stauffer. But the approach based on hierarchy 
leads to final analytical results in more complex and may be exclusive cases. 
Really, here the square approximation was taken only as an illustration. 

One has to clarify the place of the presented approach in the task to 
determine the nucleation flow. Ordinary to justify the total square approx- 
imation in the near critical zone and to use the Langer-Stauffer's approach 
one has to adopt some approximations including the smooth behavior of the 
derivative of bg near the bottom of the channel. But there is no clear evidence 
of the regular behavior of bg near the bottom. So, the approach based on the 
hierarchy is preferable. 

On the base of hierarchy one can also see many interesting and important 
facts: 

• At first we see that the quasi-unary condensation can not be described 
in terms of the square approach. A direct transformation of the for- 
mulas appeared in the Langer-Stauffer's approach does not lead to the 
formulas of the unary nucleation. This occurs because inevitably the 
square approximation has to be violated. So, we come to the impos- 
sibility of description of the quasi-unary nucleation in terms 
of the standard Stauffer's binary nucleation approach. 

• The next consequence of general results is the impossibility of situation 
of the inverse direction proposed by Zisterman-Berezhkovskii, where 
jlOj the Stauffer's approach meets difficulties. Really, now it is clear 
that valleys have to be directed to the origin, but not at the perpendic- 
ular direction as it is supposed in the consideration of Zisterman and 
Berezhkovskii. 
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One has to mention that the thermodynamics is rather formal and can 
give essential corrections to the initial variant of the theory if some other 
expressions for chemical potentials and surface energy are taken. Certainly, 
these expressions have to be the matter of discussion. But one can not deny 
the possibility to come to the situation where the square approximation is 
not suitable and one has to follow the approach suggested here. 

Here we suppose that these expressions are already given. They are some 
external information for the theory developed here. 

3 StaufFer's and Reiss' solutions 

The main goal in the investigation of the multicomponent nucleation is to 
get essential corrections in comparison with the already known approaches. 
For this purpose we shall examine the formulas of Stauffer and Reiss for the 
nucleation rate. 

3.1 Kinetic equation 

Consider the binary case. Introduce the Reiss' variables x, y as the variables 
when the free energy in the critical region has the form 

F = F,-x^ + y\ 

where Fc is the free energy in the saddle point. These variables can be 
obtained from z/i, z/2 by rotation and rescalin^. 

Instead of rotation and rescaling it is more convenient to introduce the 
separated variables directly. The variables k, ^ are the stable and unstable 
ones. One can come to 

d^F ^ dbgiO 

which is vanished in the saddle point. It means that the square form of the 
free energy in k,.^ variables looks like 

F = -A{k - + - + F, 

^ May be some part of the Lorenz transformation with an arbitrary parameter has been 
made. So, these variables aren't completely fixed. 
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without the cross term. Here A and B are some positive constants 

Then in the coordinates 
one gets 

F^Fc-x^ + f 

Now we shall seek for the similar variables obtained by the linear transfor- 
mations. 

The variables x, y can be obtained from vi, 1/2 by the linear transformation 

X = Cii(l/i - Uic) + C12{V2 - l'2c), 
y = C2l{yi - Vic) + £22(^^2 - Z^2c) 

(which isn't orthogonal) with the known coefficients 

C21 = [2(^)0] (^)c, 
C22 = [2(^)0] (^)c. 

The variables y practically coincide with x, y. The difference has an order 
of a small parameter. 

The estimates for coefficients cn, C21, C12, C22 are 

cii ~ 

C12 ~ 
C21 ~ 

C22 ~ «c 
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The estimates 

AUperp i ~ K^^^ ~ flit 

in positive powers of a big parameter k (or Utot) allows to use the Fokker- 
Planck's approximation. 

In the Fokker-Planck's approximation the kinetic equation for the distri- 
bution function n can be written in the following form 

dtn{vi, V2) = Widi[ndiF + din] + W2d2[nd2F + ^271], 

where Wi, W2 are the kinetic coefficients, i.e. the numbers of the first sort 
molecules and the second sort molecules which are absorbed by the embryo 
in the unit of time. Here 

di = d/di'i, d2 = d/dv2 

and dt = d/dt. The differentiation on the number of the molecules of the 
given sort in marked by the index near the symbol of the partial differentia- 
tion. 

Now we rewrite the kinetic equation in the variables a;, y. Note that 

di = ciA + C2idy 

82 = Cudx + C22dy 

where dx = d/dx and dy = d/dy. 

The distribution n{x,y) is proportional to the distribution n(i/i,z/2) with 
coefficient d{i'i,i'2)/d{x,y) and one has to take this difference into account 
in final calculations. In the near-critical region the coefficients of kinetic 
equation are approximately constants. 

To simplify the treatment one can use notations 

dxl — Ciidx, dx2 — Ci2dx, dyi — C2ldy, dy2 — C22dy. 

Then one can get the equation 

dtn = Wi{dxi + dyi)[n{dxi + dyi)F + {dxi + dyi)n\ 
+W2{dx2 + dy2)[n{dx2 + dy2)F + {dx2 + dy2)n] 



31 



Since the structure of terms like + dyi)F coincide with the structure of 
{dxi + dyi)n one can simply miss the last term and reconstruct it in the final 
expressions. Then 

dtu = Kidx{ndxF+dxn)+K2[dx{ndyF+dyn)+dy{ndxF+dxn)]+K3dy{ndyF+dyn) 
where 

= Wicj^ + W2CI2 
K2 = W1C11C21 + W2C12C22 
Ks = Wicl, + W2CI2 
To stress the hierarchy one can introduce the coefficients 

R^Ki, k^- — , 2- 
K2 Ki 

Then finally 

dtn{vi, U2) = R[dx[ndxF + da^n] - k~\dx[ndyF + dyu] + dylnd^F + d^n]] 

+k-^qdy\ndyF + dyn]] 
For R, k, q one can get the following expressions 

R^Wicl^ + W2cl2, 



k^ 



W,cl, + W2CI2 



WiCnC2i + W2C12C22 
^ {W^cl, + W2cl2){W^c\, + W2CI2) 

{WICUC21 + W2C12C22Y 

The last coefficient can be also written as 

C11C22 - C12C21 ,2 



? = 1 + W^W2{ 



W1C11C21 + W2C12C22' 



The value of R isn't important because it can be changed by the time 
rescaling. One can see the estimate 

k - 

which shows that A; is a small parameter. The scale of q is arbitrary, but one 
can outline situations where g — 1 -C 1. 
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The boundary conditions for the last equations are the following 



n/rf = 1 X ^ —1 — oo < y < oo, 

(8) 

n/n^ = X ^ 1 — oo<y<oo 

The plausible but not rigorous consideration corresponding to the solution 
proposed by Reiss is the following one 

• The main operator of kinetic equation is the last term in r.h.s. 

• It ensures the relaxation over the stable variable and the kinetic equa- 
tion becomes the one dimensional one. 

• The consideration of the evolution only over the unstable variable leads 
to the reduction of the kinetic equation to 

dMui, z/2) = Rd^[nd^F + d^n] 

The solution of the last equation leads to the results of Reiss. But in 
the cited paper of Reiss the hierarchy was not observed. Hence, the analysis 
there was less plausible. 

3.2 The influence on the characteristics of the process 

One needs the transformation of kinetic equation which conserves the bound- 
ary conditions, since the variables in the boundary conditions ([H]) are already 
separated. This transformation is the Lorenz' transformation. 
Introduce the Lorenz' transformation via formulas 

X + ay y + ax 



This transformation conserves the form of the free energy in the critical 
region: 

F = - ip"^ + T]"^ 
The kinetic equation is transformed to 

dMuu V2) = R{1 - a^)-^k~^[[{k - ay + a\q - l)]d^[nd^F + d^n]- 
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[{k — a){l — ka) — a{q — l)][9^[n9^F + + 

dr,[nd^F + d^n]] + [(1 - kaf + q- l]dr,[ndr,F + 9^n]]. 

Parameter a which has the absolute value less than 1 has to be chosen 
to vanish the cross term. The equation for the choice of a is the following 

{k — a){l — ka) = a{q — 1) 

Then 

dtn{vi, V2) = Ad^[nd^F + d^n\ + Cdr,[ndr,F + 

where 

^ = ^(1 - a^)-^[(k - a? + a(k - a)(l - ka)] 

C^^Jl- - ka)^ + ^^'"^^^'"^"^ ] 

k^ a 

The parameter of the Lorenz' transformation is given by 
a^l-ie + q- ^{k^ + q)^-Ak^]. 

Zrv 

After the decomposition at small k one can come to 

a = -k. (9) 

in the leading term. One can see that it is small. So, it is difficult to see the 
effect of the Stauffcr's consideration on the direction of the flow. But one 
can not directly put a = because there is a small parameter k. Expression 
for A will be 

A^R^-^. (10) 

The ratio 1/g is not small. So the correction to the Reiss' formula is essential. 
The direct substitution a = leads to 

^|a=0 — R 

which is the Reiss' result and it is not precise. 
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3.3 Some consequences for the binary nucleation 

The question to discuss here is the rate of the deviation of the Reiss' formula 
for the nucleation rate from the analogous result of Stauffer. 

In the derivation of the expression for q no suppositions about Wi and 
W2 have been made. At first the situation with the moderate ratio W1/W2 
will be discussed. 

As far as 

dui du2 di>2 
we see that the partial cancellation can take place in expression for q only in 

W1C11C21 + W2C12C22 

but not in 

C11C22 - C12C21 

So q is big enough to lead to result near the Reiss' formula A = R. This 
shows that the Reiss' formula is not so bad although it is not a true result. 

The precise coincidence of Reiss' and Stauffer's results takes place when 
q = 00, i.e. when 

W1C11C21 + W2C12C22 = 
The last relation taking into account flTT]) can be rewritten as 

Ok Ok 

OUl O'Z/2 

In the rough approximation corresponding to: 

• the capillary approximation itself, 

• the Gibbs-Duhem' equation in the capillary approximation 

• the negligible dependence of Vi on k in capillary approximation 
one can see that the last relation transforms to 

WiVi = W2V2 

where Vi is the volume per molecule in a liquid phase. This condition is the 
condition of precise applicability of the Reiss' result. It differs from condition 

Wi = W2 
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announced in paper [TU] analyzing the theory of Stauffer. 

It is clear that the last condition is wrong which opens a question of 
the formal validity of the Stauffer's derivation. Really, formally regarding 
one molecule of the first substance as several particles, one can attain the 
applicability of Reiss' result by such an artificial way. For condition WiVi = 
W2V2, this trick fails. 

Let us extract the conditions when A essentially differs from R. It can 
be only when 

g^l 

The last condition can be satisfied only when Wi <^ W2, Wi ^ W^- 
Namely, this situation occurs when there is the rapid component. The es- 
sential variation of the nudeation rate in comparison with the Reiss' formula 
is possible only under the hierarchy of the kinetic coefficients. This situation 
requires a separate analysis. 

As an illustration here we shall show the result in a square approximation 
of the free energy, although one can analytically prove that the existence of a 
rapid component throws the main nucleation flow away from the near-critical 
region and another approximations for the free energy have to be used. 

Under the hierarchy one can see the evident rapid component and for- 
mulas can be simplified. The simplification can be made also directly in the 
final formulas and the expression for A 

C21 

is proportional to W2- Then 

In the further considerations of this section this simplification is not used. 
3.4 Conclusions based on hierarchy 

In the post critical region one can assume the derivative on the unstable 
variable to be locally a constant and reduce the kinetic equation to 

dtn = R[d,{d, + h)- k-\d,{dy + 2y) + dy{d, + h)) + k-\dy{dy + 2|/)]n, 
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where h is the constant coefficient corresponding to the first derivative on the 
unstable variable and the values of i?,/;;, q are changed since the derivatives 
are taken now in the local current point. Renormalize the scale over the 
unstable variable as to put h = 1. Certainly, the hierarchy takes place after 
the renormalization. 

We are interested in the stationary solution and shall seek it in the form 

n = Q{x) exp{-{y - yo f) (12) 

with the constant mean value uq and some function Q{x). The derivative 
dQ{x)/dx can be neglected. Then 

[-kdyh + qdy{dy + 2y)]n = 

For yo one can get taking into account 

dyexp{-{y - yof) = -2{y - yo) exp(-(y - yof), 

dy{dy + 2y) exp(-(y - yof) = -4yo(y - yo) exp(-(y - yof) 
the following relation 

— = 2yo. (13) 

So the solution is obtained. 

Consider this solution. We see that the deviation of the rapid parameter 
is small also in the post critical region and the possible hidden parameter 
can not be extracted. 

Due to the slope of the free energy surface on k the minimum of the free 
energy in the cross section depends on the slope of this cross section. But 
since the slope of the free energy surface on k, is small the deviation of the 
minimum is small also. This deviation can be considered as the deviation of 
the mean value of the rapid variable and leads to the absence of the possibility 
to extract this variable in the post critical region also. 

The analogous method can be applied also for the near-critical region. In 
the near-critical region one can make the substitution 

n = P{x) exp(-(y - yo{x))^), 

where yo is now the function of x. One can determine yo according to 

{d^ - 2x)n = - J^, 
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there the r.h.s. is constant. Then 



dx'^ = —Jx + 2a;n. (14) 

The hnear character of the last equation ensures the hnear dependence of Uq 
on the unstable variable. As far as the flow is reciprocal to the halfwidht 
(along the trajectory yo) one can get the equation on the flow. The linear 
dependence of yo on x ensure the linear character of the transformation which 
is analogous to the Lorenz' transformation. 

This way of considerations can be applied to the more general situations 
without the square form of the free energy. Then the trajectory isn't the 
straight line and the solution is some approximation based on the hierarchy. 

The last question to solve is a real position of the near-critical region. 

When the deviation of the flow from the steepens descent situation is 
essential there is the danger to violate the square form of the free energy. 
The boundary conditions for kinetic equation in the critical region in reality 
have to be observed at 

n/rf = 1 -l<y<l, 

n/rf = X ^1 -1 < y <1. 
After the Lorenz' transformation 



n/rf = 1 -1<?7<1, 

n/rf = -I <ri <l. 

Rigorously speaking one has to put the equilibrium conditions at the 
line where F = Ff. — 1 which is invariant to Lorenz' transformation. But 
actually, to ensure the finite relaxation time and the constant values of kinetic 
coefficients one has to cut-off the tails and to go to the boundary of the 
reduced near-critical region. But this boundary is not invariant to Lorenz' 
transformation. 

The last definition of the boundary conditions has to be considered as the 
main one. 

But here the reduced near-critical region is stretched along one of the 
lines F = Fc where the transition occurs. The square approximation in such 
stretched region can be invalid. 
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3.5 Conclusions 

The main new results of the consideration made above are the following: 

• The hierarchy of terms in kinetic equation is shown. Earlier the hier- 
archy was observed only for halfwidths of the near critical region |25j . 

• The plausible way to derive the Reiss' formula was demonstrated. Since 
this formula is wrong, this deviation demonstrates the impossibility to 
neglect in kinetic equation all terms except the main one. 

• The moderate value of the error made by Reiss is established. Earlier 
there was a strong conviction that the error of the Reiss' approach can 
be enormous, which was illustrated by numerical examples in |B] . Now 
it is clear that the big error can be only in the cases of strong hierarchy 
between kinetic coefficients {Wi ^ W2, W2 ^ Wi) when the nucleation 
flux goes mainly far from the saddle point. 

• The simplified relations for a (see (|9])) and for the nucleation rate (see 
equation flTU]) for A) have been derived. 

• The super-critical region is studied and the expression for the distribu- 
tion function over the stable variable ( fT2l) . ( fT3l) in this region is derived. 

One can see that the precise result is rather complex. It can not be 
achieved by a simple superposition of naive solutions based on hierarchy. 
One has also to mention that even in hierarchy Wi ^ W2 the result differs 
from the naive one. 

But the main result is the absence of the really important corrections in 
comparison with a naive approach. All obtained corrections are rather small 
and mainly less than the microscopic corrections in real situations. Below, 
we shall seek essential corrections in the case of hierarchy. 

4 Nucleation rate in the situation with hier- 
archy 

The case of hierarchy certainly requires a special consideration going outside 
the local approximations in the neighborhood of a saddle point of the embryos 
free energy. 
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There are many substances for which the densities riaoi have the different 
orders of the values. For example, 

'nooH2S04 

in the everyday thermodynamic conditions. 

Assume that there are two groups of substances: the substances with a 
slow exchange and the substances with a rapid exchange. Suppose 

« W+ 

The components of the first group will be marked by the index "a" and the 
components of the second group will be marked by the index "b" . At first we 
shall consider the situation of two components and later the generalization 
will be evident. 

Here the variable v is the following one 

a 

4.1 Direction of a flow 

Extract the conditions when the flow is parallel to We construct a simple 
model which will show some estimates. 

The quantity of the embryos at the bottom with a fixed slow component 
can be estimated from above by 

Nabove = Auno eyip(-Fb) 

where Fj, is the free energy at the bottom Ai> is the effective width of the 
bottom. The normalizing factor uq in some situations of the overcoming of 
the few activation barriers can differ from the standard one. That's why we 
keep a special definition for this factor. 

The quantity of the embryos in the critical region which change the num- 
ber i/a in the unit of time is 

lA^W+AunoeM-Fb) 

The flow over the ridge is Jg- So, the necessary condition is the following 

W+AunoeM-Fb)<^Js . (16) 
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One can adopt for the following expression 
where Ft is the free energy at the top of the ridge and put Av as 



It is necessary that the transition occurs earlier than the near-critical region 
is attained. Then it is possible to put 

F, = F,-l 

at the boundary of the near-critical region. At the same boundary one can 
also put 

Ft = F, + l 

The inequality f|T6l) comes to 

« exp(2)S)-^-/' • 

Practically the same condition can be obtained by the comparison of the 
characteristic time between the transitions of the embryo along Ua which is 

and the time of the relaxation in the bottom 
4.2 The normalizing factor 

Here we shall see that there is no equilibrium distribution in the whole pre- 
critical region. 

Extract the condition when there will be the equilibrium distribution at 
the level with the fixed Va of the pre-critical region. The quasi equilibrium 
distribution has the form 

n = n'^ = nlexp{-F{ua, i^b))\ua=const ■ (20) 
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The normalizing factor Uq differs from the standard normahzing factor 
because there is an equihbrium along the band but there is no equilibrium 
between bands. 

To establish the equilibrium it is sufficient to have the intensity of the 
contact between the neighbor bands greater than the intensity of the over- 
coming over the activation barrier. So, it is necessary to determine the height 
of the activation barrier. Choose as z/^ the value of z/^e, corresponding to the 
minimum of the free energy at the band 

z/fee : miny^F{ua, I'b) = F{ua, I'b) ■ (21) 

Then the intensity of the contact can be estimated by W^uq exp(— F(z/a, z/^e). 

One can due to (fT9l) assume that the transition to the post critical region 
occurs along z/;,, i.e. inside the bancl^. Beside Vhe one can introduce u^x as 
the point inside the band where the free energy has the maximum 

max^^F{ua, v^) = F{Va, l^hx) = F^{va) ■ (22) 

Under the square approximation the transition along Uh can not occur 
because this variable is the stable one. Then i>h x can not be defined. But 
if the component Uh is supersaturated over the pure plane liquid then the 
condensation into the pure liquid is possible and z/f, ^ must exist. This shows 
that the square approximation can not be used here. 

The transition into the super critical region can occur under the arbi- 
trary Ua- But the probability of such transition is very low for all z/^ when 
exp(— Fa,(z^a)) strongly differs from exp(— F^), i.e. out of the critical region. 
But it can be greater than the intensity to come to the next band. The 
intensity of the establishing of the equilibrium (not the quasi equilibrium) 
at the next bancl§ is less than the intensity of the transition over the ridge. 
This intensity is given by 

J=Js = nlexp{-Fxiiya))W,yAxU,TT^ , (23) 

where 

W,l = W,Hu^,ua A.-. = (^^^^)-^/\.^... . (24) 

^ The value Ai/i, depends on weakly. 

^ The intensity of transition to the next band. 
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There is no need to establish the equihbrium along the whole band with 
the small z/^. The value of z/^^. for small z/^ can be very big, the barriers of the 
nucleation can be very high, but it is necessary to have the equilibrium only 
near the bottom, i.e. at near to Ube- The establishing of the equilibrium 
along the whole pre-critical region of the band is necessary only for the bands 
where the intensity of the transition to the post critical region is essential 
(comparable with the intensity of the transition between the bands). Ac- 
cording to the previous considerations there is the quasi equilibrium along 
such bands. 

Introduce the number of embryos in the band 

N{i/a) = nlexp{-F{iya,n e))Ki^b , (25) 

where 

^e^b = ^'^V{-F{iya, Ub) + F{Ua, Ube)) (26) 

has a sense of characteristic width. The last formula in a continuous limit 
can be transformed to 

A^Ub^ / ex.p{-F{ua,Ub) + F{iya:iybe))dn ■ (27) 
Jo 

At the ends of the interval of integration the equilibrium distribution 
can be violated but there the subintegral function goes to zero. As far as 
exp(— F) as function of Ub is rather sharp near the maximum then the number 
of the embryos going from the band with z/^ to the band with z/^ — 1 can be 
approximated by W~{ua,i'be)N{i'a). The number of the forward transitions 
is W^iva ~ 1) ^be)N{^a ~ !)• Then one can write the balance equation at the 
band 
dN 

= W^i^a - 1, l^be)N{Ua - 1) + {u, + 1, Ube)N{u, + 1) - 

(28) 

For J{ua) one can get 

J(z.„)^Ar^^-^^P^-^^^°'^^-^--^-^^"^^ . (29) 

^eJ^b^xJ^b 

One can see that the absence of the equilibrium distribution in the whole 
pre-critical region is the characteristic feature of the transition far from the 
saddle point. 
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4.3 Valley zone and ridge zone 

For every Ua in the pre-critical region there will be i/^e- The curve Ubei^a) will 
be called the valley in i/^, vi, plane. 

For every Ua in the region under consideration there will be i/bx- The 
curve Vhxi^o) will be called the ridge in z/^, plane. 

Since there is a slope of the ridge and the valley in t/^ direction it is 
necessary to specify the set of variables. 

In the set of variables k, ^ the channel of nucleation is the straight analog 
of a valley. But the channel of nucleation does not coincide with the the 
valley in Va-, fh plane. 

The line analogous to the ridge, i.e. the ridge in ^ plane will be the 
separation line defined as 



dFiK, 





< 



The values at the channel of nucleation here will be marked by the sub- 
script h and at the separation line the values will be marked by the subscript 
s. 

We see that effectively the flow is directed along uj,. The problem to 
get J(t'a) is purely a one dimensional problem. So, in the band = const 
there exists the valley Vb ^ e zone and the ridge I'b ^ x zone. Precise 
definitions are the following 

• The ridge zone in i/b scale is determined by conditions 

Certainly, F{i'a,i^b) < F{va-i^bx)- This zone has to be near the given 
ridge. 

• The valley zone in Vb scale is determined by conditions 

Certainly, F{ya-i^b) ^ F{i'a,J^be)- This zone has to be near the given 
valley. 
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To find the value of ttie flow J(z/a) one has to solve kinetic equation in the 
ridge zone. To find the normalizing factor like it was done in heterogeneous 
nucleation it is necessary to consider the valley zone and to solve kinetic 
equation in this region. 

The problem under consideration is the influence of the surface excesses on 
the forms of the free energy in the ridge zone and the valley zone. Fortunately 
some simplifying properties will be established below which help to escape 
from the explicit inclusion of surface excesses in the kinetic equation. 

For the ridge zone these properties are the following 

• Define by the subscript the values without surface excesses 

• In the ridge zone for arbitrary s corresponding to the ridge zone 

For the valley zone these properties are the following 

• In the valley zone for arbitrary s corresponding to the valley zone 

One can analogously define the channel zone and the separation zone. 

• The separation zone is determined by conditions 

F{k,^)>F{k,Q-1 

The value of k is fixed here. Certainly, F{k,,^) < F{k,,^s)- The sepa- 
ration zone has to be near the given separation line. 

• The channel zone is determined by conditions 

FiK,0<FiK,^n) + l 

The value of k. is fixed here. Certainly, F{k,,^) > F^k.,^^). The channel 
zone has to be near the given channel line. 

One can analytically prove the following properties for the separation 
zone 
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• In the separation zone for arbitrary s corresponding to the separation 



One can analytically prove the following properties for the channel zone 
• In the channel zone for arbitrary s corresponding to the channel zone 



The method of a proof of all these properties is quite analogous to the 
already presented for the near-critical region. These properties allow to solve 
kinetic equations in these regions by some shift renormalizations and solu- 
tions in the absence of the the surface excesses. 

4.4 Discrete case 

Consider the stationary solution. The last equations form the system of 
algebraic equations. Note that the sufficient equations arc those where W^N 
has the order of J. The equations with W^N ^ J can be taken into account 
by the boundary condition n = n'^ = for Va which is less than some Vamin-, 
where J begins to be comparable with W^N . More precisely this question 
will be discussed later. 

Formally one has to put this condition at ^ i^amin- Then one has to 
solve equations and to see where the condition n ^ rf will be violated. It is 
very easy to do having calculated J{va) on the base of rf to get 



zone 



F(«;, is + s)- F{k, Q ~ F{k, U + s) - F{k, U) 



F{k, ^h + s)- F{k, Ch) ~ F{k, ^0 + s)- F{k, Cho) 




or having expelled the unphysical region 




This will give the necessary estimate. 

In the region where W^N <^ J the solution is rather simple 



(30) 
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This condition will be seen automatically at some i>a and since n/n^ is a 
decreasing function of Ua it will take place later. So, one has to investigate 
only few equations of the type 

W^iUa - 1, iybe)NiUa - 1) + + 1, Z/fee)iV(z/„ + 1) 

(31) 

The total flow is defined as 

Jint= E J^^a) , (32) 

where I'amax marks the upper boundary of the equations sufficient for the 
consideration. 

In the limit when there is only one sufficient equatioE0 

Jint = W^ noexp{-F{ua,i^b))dub =W^Ntot{i^a) (33) 
Jo 

where the total number of droplets at Ua is 

Ntoti^a) = noexp{-F{iya,iyb))diyb 
Jo 

The discrete situation is the most frequent one. But namely this situation 
has not been considered earlier. 



4.5 Differential model 

Consider the opposite situation when among fl3T|) there are so many equations 
that it is difficult to solve the algebraic equations. Then it is reasonable to 
come to the differential form. The condition of the validity of the differential 
form coincides with the condition of the big number of the essential equations. 
Then 



(34) 



rj2 

{iW^iUa,Ube)+W-{Ua,Ube))N{Ua)} 



^ Having attained i^a all embryos come automatically to the super critical region. Then 
it is possible to write the expression for the transition on i^b. 
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With account of (12^ one can get 



^ {{Wg^u^, - exp(^^%^)))iV(z.J} + (35) 



One can note that 

• The hierarchy of the halfwidths of the near-critical region shows that 
the quasi-unary nucleation in the square approximation in the neigh- 
borhood of the saddle point is impossible. 

So, the change of approximation to a linear one is absolutely necessary. This 
conclusion is very essential for further consideration. 
One can use the following approximations 

J = Joexp(cy) , (36) 

y = Ua-UaO , (37) 
_ dF{Ua,Ube) : dF{Ua,Ub^) 

It means that the linear approximation for F{va-, i^bx) — F{va-, Vbe) is adopted. 
The supposition made in this paper radically changes from the supposition 
of Trinkaus. This difference will be discussed in a special part of this paper. 
One has to note that 



differs from 



and 



dFjug, z/fc) 

dF{Ua, Vbx) 
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differs from 

When we use '^'^ ^9^'^'" ^ imply the differentiation along the bottom of a 
valley. When we use ^^^gl'a'''^ 

imply the differentiation along the top of 

a ridge. 

Then one can get 

IeMcy)N = - 6)— + W^il + 6)^ , (40) 

where 

W+ = Wt{ua,ybe) , (42) 

e = exp(^-^) . (43) 

It is supposed that e depends on Va rather weakly. We suppose that e is 
locally a constant value. This supposition is many times weaker than the 
previous approximation. 

Since is small the value of e is close to 1 and 1 — e is very small. 

Then the value 1 + e is close to 2. Then the relative deviation of ^^^a'"'^'"''* 

ava 

have no importance. 
Then one can get 

x^cy , N exp(x) + A-— + B— = (44) 

dx'^ dx 

with the known values of ^4, B. 

After the transition to = exp(x) one can get 

A^'^N" + (A + B)i)N' + 7pN^0 (45) 

with the known solution 

Ar = ^-^/(^^)Z«/4-^^V2) ^ (46) 

where Zj is the cyhnder function. One has to choose the solution vanishing 
at 00. 

The known value of allows to determine the total intensity of the 
embryo formation and the integral can be taken analytically. 
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4.6 Applicability of solution 

Our solution corresponds to the solution derived by H. Trinkaus in [9]. But 
this correspondence is only a formal one. Recall the derivation by Trinkaus 
in [9]. Trinkaus proposed the linearization of the free energy F {G in terms 
of Trinkaus) around n2 (this value is analogous to i^ao)- 

Now we shall analyze the possibility of linearization of F in the vicinity 
of Uao- This linearization can be considered in the global sense and in the 
local sense when linearization is done over one coordinate while the other 
coordinate determines the values of coefficients in this linearization. 

Linearization in the global sense can not exist because the second deriva- 
tive at the ridge and the second derivative at the valley must have different 
values. Only then the value of 

AF{iya) = F{Ua, Z/fe^) - F{Ua, Vbe) 

will be a real activation barrier. The exponent of the last value is the leading 
term in the expression for the flow. 

Linearization in the local sense can not be valid also. It is absolutely clear 
that the linearization over z/;, can not be made because it is necessary to have 
a valley and a ridge for F as a function of z/;,. So, it can not be linearized. 
Another possibility is to fulfill linearization over z/^ while coefficients depend 
on z/{,. The last possibility is the most preferable one. 

The careful analysis of the last possibility shows the impossibility of lin- 
earization. Really, since the ridge in z/^, z/f, scale is relatively close to the 
ridge in k, ^ scale one can see that the behavior of F as a function of Ua at 
z/fe slightly greater than z/fej. is the following one: At first F increases until the 
ridge in z/q, will be attained. Later with increase of z/fe the value of F will 
decrease. This behavior is the direct consequence of the slope of the channels 
of nucleation in z/^, z/f, plane. So, the linearization is impossible. 

The only possible variables, in which the approximate local linearization 
is valid are variables E,- One can see that there F can be linearized far 
from the critical point 

dF{K,0/dK = 

at every C,- The linearization is made only along n. But these variables have 
not been even mentioned in [9]. 

It has been already analytically shown that we are far from the critical 
point. Namely this allows the linearization in a local sense along k. 
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The critical point which is the nearest to the origin of coordinates is 
situated in the channel in z/^, z/5 picture. This is the real saddle point. But 
since we are far from the main saddle point it means that we are far from 
every critical point. 

Now we shall see that the linearization of the free energy in k, ^ variables 
is possible. Really, 

The second derivative is 

The size of characteristic region in which the linearization is necessary 
can be estimated as 

Ak = (6,(0 - 2«:-V3/3)-i 
So, the necessary condition is 

|(-6,(0 + 2«:-^/V3)-22«:-^/V9| < 1 

Since we are far from the critical point one can neglect the compensation in 
{-bgiO + 2«;"^/V3) and get 

\{2K~^/y3)-^2K-^/y9\ < 1 

or 

« 1 

The last inequality is evident. 

The last property is important for our needs. We are interested in the 
linearization of the free energy of the ridge and of the valley. Really, the 
particular case of the last derivation is the possibility of linearization of F 
along the ridge and the valley in k, ^ scale, i.e. along the channel and along 
the separation line. 

The last step is to go from k, $, picture to Uf, picture. We see that the 
slope of the valley and the ridge in k, ^ picture along k is very small. Since 
the slope is proportional to | ^^q*^'^^ \ it can be seen from 

l^^^l = I - b,{0 + 2K-'/y3\ ^ 2/.-V3/3 « 1 



51 



So, the characteristic distance where the height of the vaUey, the height of 
the ridge and, thus, the height of the activation barrier (in fact it can be 
proven that there is no compensation) undergo the variation of one thermal 
unit is 

One can see that Di <^ k and it means that the relative size of the transition 
region has to be small. 

This slope has to compared with the characteristic halfwidth along i/h or 
the characteristic size Di has to be compared with the half-width along ^ 
multiplied on k. We have 

We see that 

D2<^Di 

The slope at the boundary of halfwidth is 

and it is rather essential. 

We introduce the distance D3 where the slope 



has the order of the slope of the ridge dF/dn ~ k, i.e. k, Then we 
get 

D3 - «;2/3 

We see that the order of D3 is the same as the order of Di and it is relatively 
small 

It means that the deviation of the separation line in k, ^ scale from the ridge 
in scale is relatively small. 

Since F^, Fg allow linearization as functions of k or of Ua we come to a 
conclusion that the linearization of F^, Fr, AF (this value is a function of 
one variable) as a function of k or of i/^ is quite possible. 



52 



4.7 Simplified solution 

Since dF{va, Vbe) /dva ^ 1 one can put e = 1. Then 5 = and one come to 
the universal solution 

N ~ Zo(4x^'^') (47) 



This is the universal function Zq of the variable 

2 

exp(ca;/2) 



'A 

Finally we get a universal solution 



4.8 Discussion 

The multidimensional case is quite analogous to the two-dimensional one. In 
the multidimensional nucleation one has to consider some channel of nucle- 
ation. One has to extract the set of fast variables {v^} and the set of slow 
variables {uq}- 

For the set {ua = fixed} one can establish J{ua} by the consideration of 
the evolution in the set {i^b}- It can be done by the standard methods from 
the previous sections. 

After the calculation of J{i,a} one can define the direction. It will be the 
quasi- integral on Ua- This defines the first coordinate. The second coordinate 
is the direction of the bottom of the valley in the cross section {uh = const}. 
The further consideration is absolutely analogous. 

The new results formulated above are the following: 

• In the paper of Trinkaus [9] only the differential case was considered. 
The discrete case was not considered there. Really, the height of the 
pseudo-activation barrier can change rather rapidly with increase of Ua. 
This leads to the preference of discrete model. 

As for the half-widths of the bottom of the channel and of the top 
of the ridge in calculation of J there are inequalities which guarantee 
the possibility of the differential description. Really, these half-widths 
increase like n^^'^ (see the standard estimates for the half-widths along 
the stable variables). But if even these variables will be not so big 
nothing will be changed because they variate slowly in comparison with 
the exponent of the height of the pseudo-activation barrier. So, the 
mathematical structure of the balance equation will be the same. 
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• Here the surface hmited growth is considered while in [H] the diffusion 
hmited growth was used. It seems that because the transition occurs 
earher than the saddle point will be attained the embryos are small 
enough and the surface limited growth is preferable. 

• It is shown that the absence of the equilibrium distribution in the pre- 
critical region is the driving force of the transition far from the saddle 
point. This fact stresses once more the importance of the formulation 
of the boundary conditions and outlines the paper [8] where the bound- 
ary conditions were used for the situation without hierarchy of kinetic 
coefficients. 

• The hierarchy of the halfwidths of the near-critical region (more ac- 
curate the near-saddle region) shows that the quasi-unary nucleation 
in the square approximation in the neighborhood of the saddle point 
is impossible. So, the change of approximation to a linear one is ab- 
solutely necessary. Moreover, it is impossible to see the transition of 
the Stauffer's solution to Trinkaus' one on the analytic level of explicit 
formulas. 

Beside the mentioned disadvantages of the differential approach one can 
mention the disadvantage connected with the position of the basic point u* 
for decompositions of the height of the ridge and depth of the valley. An 
ordinary chosen point for such decompositions is 

W+ = J{va)/N{v,) (48) 

The presence of this point awakes the idea of the Genuine Saddle Point |TT] . 
It is reasonable to put the point of decomposition at 

n{K) = ne.K)/'^ (49) 

The shift between i/* determined by (HSj) and will be called "the soft 
shift". 

The greater is |c|~^, the greater is the soft shift. But the applicability of 
differential approach requires 

|c| < 1 

The last parameter ordinary comes from two decompositions: one of the 
height of the ridg^ 

*Take a cross section {i/b = const}. 
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with parameter and another of the depth of the vallejo 

with parameter k^. 
Ordinary 

ke>0 

(the opposite sign means that the saddle point is aheady behind) 

kr <0 

(the opposite sim means that energetically it was more profitable to cross 
the ridge earlieio). Then in 

J =i Joexp(-A;,.(z/a - Uao) + kei^a " l^ao)) = JoeXp(c(z/„ - Uao)) 

parameters kg and kr can not be compensated. Ordinary both linear approx- 
imations are necessary. 

Then the condition |c| <^ 1 leads to 

\kr\ < 1 
\ke\ < 1 

Under the last two inequalities one can see that N becomes many times 
less than the equilibrium value A^*^'' much earlier than z/^ = Uao and the 
transition is actually over. So, the point of decompositions has to shifted. 

The shift of decompositions has to lead to the basic point situated at the 
position characteristic for the relatively intensive flow. One of the possible 
recipes is to choose the point z/* of decomposition according to 

Jo 

The last condition can be approximately rewritten as 

\K\ + m N 

^Take a cross section {j/f, — const}. 

^°Then the cross of the ridge can not disturb the equihbrium distribution. So, the flow 
is known. 



55 



One can start instead of fa = from infinity and get a similar estimate. Also 
it is reasonable to consider 



|A;^| + |A;e| 2N 

as the point for decompositions. 

Here naturally appears the length A of the region where the transition 
occurs. It can be estimated as 



Jh-p ~\~ }Zp 



So, the soft shift can be greater than this region. 

We continue to consider the problems of the differential approach. 

• Another problem is the smallness of | /cg | , \kr\. Because of the monotonous 
character of derivatives of the free energy along channels and ridges it 
can be attained only near the saddle point. But here the square approx- 
imation has to be used and the Stauffer's solution will be the answer. 

Certainly, if the value of is extremely big one can observe small 
values of derivatives rather far from the saddle point. But, although 
even here the discrete approach is preferable as it will be shown later. 

Now the simplified approximate method for continuous case will be pre- 
sented. In equation 



dx^ ^ dx Az/,Az/, W+ 



one can put to zero because of the smallness of \ke\. 

Also because of the smallness of |A;e|, \kr \ one can put very approximately 

exp(F, - Fe) Ty+ 
Ai/eAi/, W+ 

to some constant (let it be /q)- Then 

dx"^ ° 
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Solution of the last equation is evident 

N = Aex.p{-^/I^x) +Bexp{^x) 

The requirement A?" — > at x — > oo leads to 

N = Aex:p{-^JTox) (50) 

But this solution has a bad behavior at x — > — oo. So, in this region one 
has to use another approach. At a; — > — oo the flow is very small and N is 
approximately equal to the equilibrium value Neq. Then 

m^k'^^N exp(fi-Fe) iy+ 
^dx i^v.i^v. 



Then approximately 



iV = iVeg - / Jdx 



or 

TV = 7V,„ - / TV, """^ V"'".' '~'''"'~' dx 



_ f eM-{\ke\ + \kr\)x) 



With the evident approximation for the equilibrium value Neq'. 

Neq = 7V*exp(-|A;e|x) 
with parameter — Neq{x — 0) one can get 



*J AUrAUe W+ 

Since one can approximately take Ai/^Ai/g as a constant value there are no 
problems with integration. So, 

^^^^^_^^l-exp(-|A..|x)Ty,t 



IkJAurAu^ W+ 



These two solutions (let it be A^i and A^2) have to be stuck together at 
the point where 
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Another method can be formulated if we notice that flSUI) is vahd namely 
locally because it was derived with a supposition Jq = const. So, we have to 
go to the local form by differentiation of fl50l) which gives 



This equation can be integrated with arbitrary Jq which leads to 



is taken, one can come to 

~ exp( — — exp(cx)) 
c 

with parameter Iqq. Certainly, parameters /qo and c can be considered here 
as the fitting parameters. 

The functional form announced above resembles 0-function with a soft 
transition from 1 to 0. We shall call it as a soft 0-function and denote it by 

S{x) = exp(— exp(x)) 

This function can be used brick in an ansatz 



which can be very effectively used as an approximate solution in all situations 
considered below in this paper. 

5 Interaction of valleys 

5.1 Coordinates of valley in the iya,i^b coordinate sys- 
tem 

The coordinate of the valley is given by the condition 



A^ ~ exp( 




When the evident known functional form 



Jo ~ exp{cx) 



Q = "^AiSiaiix - Xi)) 
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The straight differentiation of the free energy gives 

dba dbh 

where S is the surface square of the embryo. In simplest approximation it 
can be written as 

S = {vaUa + vmfl^ 
The standard Gibbs-Duhem's equation looks hke 

Vadha + Vbdhb = 

and leads to 

dha dbb dC. „ 

This brings the condition for the valley coordinate to 

dvb d^ dvb 3 ^ ^ dvb ^ d^b ^ 

But due to the surface enrichment the concentration differs from 



l^a + T^b 

and has to be 

^ _ l^a - SQa 

T^a - SQa + l^b - S Qb 

Then the Giibs-Duhem's equation looks like 

Sd^ + Vadha + Vbdhb = 

and leads to 

+ \—v — + —V —1 = 
d^ dvb dvb dvb 

and 

9F 2 i/2r , ^^'^ 9^ , 
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The careful analysis of the generalization of the Gibbs-Duhem's equation 
for the embryos shows that the terms 

have to vanish together with 

dub dC, " dub dC, dub 
Really, the Kelvin's relation in the saddle point requires that 

ba _bb 

Va Vb 

The direct calculation with a non zero value of the last terms gives 

2 n-l /2 ba 



(51) 



7(367r)V2f5-i/: 



U,, _|_ dva,. I dvt, d^i ^ \i) -i- ^^"u -I- ^^fc/A ^ 

[Vb + -dl^ag^^ + ^ [Va + af'^'^^^ + 'W'^^d^^ 

and one can come to flSTj) only if these terms vanish. 

Generally speaking the Gibbs-Duhem's equation has to written in the 
form 

y^{differentials of all intensive variables)*{cor responding intensive variables) 
Particularly 

5^7 + Vadba + Ubdbb + ^adVa + I'hdVb = 

Then 

dj ^ 2 y^ dva dvb dba dbb 

and 

^-^ 2 ^2 

_ = / Vb - bb 

OPb 3 

One can see that the concentration of valley satisfies the condition 

and it is not a constant value. Moreover, it is evident that valleys in the 
^ai system of coordinates do not coincide with channels in /t, ,^ system of 
coordinates. They coincide only in saddle points. The valleys in Vai system 
can appear and disappear, their position in the absence of hierarchy of kinetic 
coefficients means nothing. 
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5.2 Asymptotics at — ^ oo, z^^ — fixed 

The necessary condition of applicability of solution of Trinkaus is the limit 

1/}, 00, i/a — fixed F —00 

The explicit calculation gives 

i/h -^00, fixed F -bbCb 

So, it is necessary that Cb > 0. But the last condition is not a necessary 
condition for nuclcation in a gas mixture. The necessary condition is the 
existence of concentration ^ for which the function ba^ + bb{l — ^) is negative. 
So, there exists a situation when there is no behavior necessary for application 
of the Trinkaus' solution. 

When bb > one can see the asymptotic wing with a negative slope. It 
will be called simply as "wing". 

Otherwise there can be a situation when even with 65 > nucleation 
can go from one valley to another (may be more deep) valley and further no 
transition to 1/5 —> 00, I'a = fixed will take place because of the height of a 
new (further) ridge. 

In the case of purely supersaturated vapor of components the wings have 
to be included in the general picture of relief of the free energy. 

5.3 Two valleys. Kinetic equation 

The previous consideration shows that the most ordinary situation is the 
jump of embryo from one valley to the neighbor one. In one valley (let it be 
called as the "source valley" and marked by the subscript — ) the embryos are 
in the pre-critical region (i.e. k < Kc) and in the other valley (let it be called 
as the " destination valley" and marked by the subscript +) the embryos are 
in the post-critical region (i.e. k > Kc)- The transitions take place along 
lines Ua ~ const. Since the increase of i/b leads to the increase of k it is quite 
possible. 

The values of Uf, at the ridge will be marked as Ubr- The values of i/f, at 
the bottom of the source valley will be marked as ly^e- and the values of i/b at 
the bottom of the destination valley will be marked as i'be+- All these values 
are taken in the i/a, coordinate system. 
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Kinetic equations are rather transparent and look like 
dN 

= W+{Ua - 1, l^be-)N4Ua - 1) - W+{Ua, l^be-)N.{Ua) 



dt 

(52) 

+W7.(^a + 1> l^be-)N4Ua + 1) - i/fte-) A^-(i^a) " ^-(i^a) + ^+(i^a) 

^ = W^a'd(^a - 1, ^6e+)A^+(^a - 1) - W^<i(^a, ^;.e+) A^+(^a) 

(53) 

+W^{^a + 1, i^6e+)A^+(j^a + 1) - W^{Ua: iybe+)N+{Ua) - J+{^a) + J-{^a) 

Here and are the numbers of embryos with given Va in a valley (in 
^a, system of coordinates), and VF_ are direct and inverse absorption 
coefficients, J_ is the flow from the source valley to the destination valley, 
J_ is the flow from the destination to the source valley (the inverse flow). 
We shall investigate the stationary solution. 

One has to take into account that and W~ are functions of Uf,. They 
are taken at Ub equal to the values at the bottom of valley. This can be done 
because of the relative narrowness of valleys which goes from representation 
in ^ variables. 

5.4 Two valleys. Direct and inverse flows 

The values of flows J_ and J+ are given by the standard formulas 

J ^^exp(-fi + F._) 

ArUbK-J^b 

J - M exp(-F, + Fe+) ^+ 

Arl^b^e+l^b 

Here Fr is a free energy of the embryo at the ridge (in Ua, Vb coordinates), 
Fe_ is the free energy of the bottom of the source valley in Ua-i ^b coordinate 
system, Fe+ is the free energy of the bottom of the destination valley in i/^, 
Vb coordinate system. 

The value of A^Ub is given by 

Ar^b = CXp{-Fr + F{Ua, Ub)) 
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Here 1/5^1 a-nd i^6r2 are chosen as roots of equation 

F{Va,Vb) = (2F, + Fe+ + Fe_)/4 

closest to viyr and 

^brl < ^hr < ^br2 

The value of Ae-Vb is given by 

Ae-Ub = exp(Fe_ - F(l/„, Ub)) 

Here I'be-i and i'be-2 are chosen as roots of equation 

F{iya,iyb)^{Fr + F,_)/2 

closest to Ube- and 

The value of A^j^Ub is given by 

Here Vbe+i and fbeH-2 are chosen as roots of equation 

closest to Vbe- and 

In continuous approximation one can get the following equations 



Arfb = / exp(-Fr + F{Pa, i^b))di^b 

r^bB-2 

Ae-Ub = / exp(Fe_ - F{Ua, l'b))dl'b 

r^be+2 

Ae+Ub = / exp(Fe+ - F{Ua, Ub))dl'b 

One can prove that in frames of inequalities lying in the base of capillary 
approximation the continuous approximation is valid. One can also prove 
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that in the absence of pecuharities in behavior of the free energy it is possible 
in the capillary approximation to use the square approximation with infinite 
limits for calculation of the mentioned values. This gives 



A . ^-l/2 

A^z/fo = V7r(-- W=vw) ' 

Ae_i/;, - V7r(- W=v^,_) 



One can rewrite equations for J_ and J+ as following 

where 7+ and /_ are independent on N+, A^_. 

Already now one can fulfill the qualitative analysis of the kinetic equa- 
tions. 



5.5 Qualitative analysis of the kinetic equations 

Consider the region of Ua where W^{iya, Vbe-) ~ It is easy to see that at 
Ua corresponding to the possible transition from one valley to another 

(otherwise the saddle point in the source valley is already over) 

(otherwise it will be necessary to overcome the saddle point in the destination 
valley and it will cause the establishing of the equilibrium distribution until 
the height of the saddle point; moreover there is a straight way without 
barriers to the origin of coordinates). 
It means that 

F{Va, l^be+) > F{Pa + 1, l^be+) 
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Moreover one can see that 



(otherwise it is more profitable to overcome the ridge earher at smaller i/g)- 
Practically in the main order 

^— = exp(+Fe_ - Fe+) 
-'+ 

The ratio governs the evolution of the process. One can see two 

characteristic situations here 

• Situation 

/- >/+ 

Here one can see that the solution of the previous section can be directly 
applied. The destination valley do not affect the distribution in the 
source valley. So, one can put J+ = and split the system of equations. 
Only the first equation is essential and solution is really the solution in 
the situation discussed above. 

• Situation 

/- </+ 

This situation has no analogs and has to be considered separately. 
5.6 Situation 

/-«;/+ 

One can approximately put 
This is taken only for simplicity. 

Approximately, the condition of the beginning of the jump of embryos, 
which changes is the following 



I->W+{Ua,l^be-) 
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Then 

Then the second equation of the system becomes the following 

J_ = J+ 

and we have locally in a rough approximation 
dN 

= W^{Va - 1, V,,_)N_{Va - 1) - W^{Va, V,,_)N_{Va) 



dt 



(54) 

+W-{Ua + 1, Ube-)N4Ua + 1) - W- {Ua, i/foe-) A^- (i^a) 



Then 

IT— = ir^ exp(-Fe_ + Fe+) 
and approximately 

^ = exp(Fe_ - Fe+) 

The point where 

will be marked as Ua — yo- When Ua increases one has 

I+:$>W+{Ua,l^be-) 

This ensures the quasi-equilibrium and actually the common valley. 
Later one attains yi where 

Fe_(yi) = 

For Ua > Hi one has 

/_ > /+ 
J_ = J+ 
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It will be until 1/2 defined by condition 



/+(7/2) = W^il^a = 2/2, l^be-) 

(also the soft shift has to be added). Later all remaining embryos from 
the source valley go into the destination valley. But their total quantity 
is already rather small. So, we need not to consider this process in 
details. 

The main conclusion results in the appearance of the common valley 
with a new free energy Fq. Here there is no connection with the absence 
of excesses. This free energy can not be defined separately from the 
width of the equilibrium distribution, only the ratio 

exp(-Fo)/Aeoz^6 

can be determined. But namely this ratio is the equilibrium distribu- 
tion and in the expression for the nucleation rate. 

The last ratio can be determined from 

exp(-Fo) ^ exp(-Fe-) _^ exp(-Fe+) 

Aeot'b Ae-l^b Ae+I^b 

Very approximately one can say that 

exp(-Fo) ^ exp(-Fe_) 

when Fe_ < Fe+ and 

exp(-Fo) exp(-Fe+) 



when Fe_ > Fe+. 

5.7 Intermediate situation 

Intermediate situation is very rare because it can take place only under 
the simultaneous realization of two equations 
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and 

(also the soft shift has to taken into account). But this case in the 
only one where the interaction of valleys and the exhaustion of the 
equilibrium distribution play simultaneously. 

Solution of this situation is rather simple - it is necessary to solve the 
system of several algebraic equations. At small 

< yo 

where 

one has to use the boundary condition 

N- = N^eq ~ exp(-F(i/„, i/be-)/Ae-i^6 

< iV+eg ~ exp(-F(z/„, Z/be+)/Ae+Z/b 



At big 
where 



J^b > yo 



W+{iya,iyhe-) < I- 
one has to use another boundary condition 

/+ = 

if it will be necessary. So, the task is to solve several simple alge- 
braic equations. Certainly, the discrete approach is preferable in the 
computation. 

To come to the continuous approximation one has to change the finite 
differences for derivatives which approximately leads to the following kinetic 
equations 
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ON. 



dt 



+ 



+ 



+ k 



dN. 



a 



+ 



Here 



(1 - exp{dF{uaUhe-)/dua)) 



dF{UaUhe-)/dUa 



exp{dF{l'al'be+) / dl^a) 



dF{UaUbe+)/du, 



a 



Continuous approximation can not be widely spread but can be applied 
only in rather specific situations. The reasons are similar to those described 
in analysis of the Trinkaus' solution. The proximity of dF^/dvi, and both 
dFe+fdi/fj and dFe-fduf, to zero means the proximity to the saddle point 
where the linear approximation fails. 

The simple approximate method is the iteration one - the values J_ and 
J+ are calculated on the base of previous approximations and they are treated 
as known functions. Initial approximations are following: 

• when the source valley are deeper than the destination one, then there 
is the quasi-equilibrium. 

• when the destination valley are deeper than the source one, then there 
is the Trinkaus' solution or the corresponding simplified solution. 

This method is very effective and leads to a rather precise solution after 
one or two iterative steps. 

It is necessary to stress here the effectiveness of the method based on the 
ansatz with the soft Heavisaid's functions. 

The main result of this section which was the goal of the whole publication 
is the radical change of the nucleation rate. The main goal is achieved - the 
change of the nucleation rate in the orders of magnitudes is shown. One 
can also see that the rate of nucleation does not depend on the free energy 
in saddle point but on the mutual position of valleys and ridges and their 
relative heights. Certainly, the problem to find the nucleation rate includes 
now the determination of many characteristics and is more complex than 
in the theories suggesting the recipes based on the value of the free energy 
in one point. The theory presented here has to be used in order to get 
the true value of the free energy. Now the problem is transformed in the 
thermodynamic area - it is necessary to find the free energy of the embryo 
formed in the mixture of vapors. This problem is complex enough to continue 
investigations of the binary and multicomponent nucleation. 
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6 Paths of transition 



Now one can return to the general situation to see how the real transition 
from the pre-critical region to the post-critical region will occur. 

The problem is to see where the real change of the channels will take 
place. This problem will be solved here. So, here the analysis will be mainly 
qualitative. All details of transition between channels will be a subject of a 
separate analysis. 

6.1 Approximate position of the valley 

To get the approximate position of valley and the ridge one can act without 
surface excesses. 

Consider the channel in coordinates i^a,^- Then 

where p{^) is a known function and 

The coordinate of the valley is given by condition 

dF ^ 

or 

At the saddle point 

and the saddle point of valley coincides with the saddle point of the channel 
line, since 

B'iO = 

Asymptotically at — > oo one can get 

-B'iO'fiO - BiO'p'iO = 
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One can see that the function p is rather smooth while B is rather sharp. 
This condition is a definition of a "clear channel". Then one can neglect 
B{i)p\C) ill comparison with B'[^)p{^). This leads to 

-5'(0p(0 = 

and because of p 7^ the last equation coincides with the coordinate of the 
channel. So, we see that the valley is near the channel line for every Va- 

To see the behavior at moderate fa near the critical values one can note 
that p' attains a moderate value. Really 

p = ~ 1 

(we choose the space scale to have the volume for a molecule in a liquid phase 
the order of 1). Then one has to take into account that 

= 2/{3max B{^)) > 1 

requires max B <^ 1 Then the term B{^)p'{^)i'a has a small parameter. The 
term p~^/^(^)z/~^/^p'(,^) has the same order as B{^)p'{^)h'a and, thus, is small. 
This reduces the coordinate of a value to coordinate of a channel. 

The same analysis can be done for every ridge. The general approximate 
conclusion is that every separation line corresponds to the ridge and their 
coordinates are similar. 

At fa — > and k — > the leading term is 

which means that the valley does not exist. So, the valley can not directly 
start at Ua = in continuous approximation. Fortunately, ordinary this effect 
takes place at less than 1. 

All above considerations are very approximate and they are used only to 
see that qualitatively nothing is changes when we consider valleys instead of 
channels. 

Approximately speaking every channel corresponds to one valley, their 
coordinates are rather similar. 

Precisely speaking one can see many specific peculiarities, for example, 
the appearance of valleys without corresponding channels. But the proba- 
bility of such peculiarities is very low. As a rule these valleys are not deep 
enough and can be treated as negligible ones. 
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6.2 Transition zones 



Consider the pair of valleys. 

Every valley (index v) can be considered as a source valley (s). Every 
valley can be considered as a destination valley (d). One can imagine many 
pairs of source and destination valleys. Every pair has to be investigated. 

At first we consider the situation when the channels are neighbor ones. 
The ridge is the maximum of F at the band i/a — const between the con- 
centration of a source valley and the concentration of a destination 
valley. 

We define Fy as the free energy at the valley, F^ the free energy at the 
ridge. 

Now we shall make use from the approximate functional form for Fj., F^ 
estabhshed above 

F — constivl^^ — const2i'a 
One can see the following facts 

• Every valley has only one critical Vavc point determined by 

dFy/dua = 

• One can define the pre-critical region of the valley where dF^/dva > 
and post-critical region of valley where dF^/dva < 0. There is only one 
pre-critical region with a size Ua < i^avc and a post-critical region where 

^ tave- 
rn Every ridge has only one critical Uavc point determined by 

dFr/dua = 

• One can define the pre-critical region where dF^/dva > and post- 
critical region of ridge where dFr/ dva < 0. There is only one pre-critical 
region with a size Ua < Varc and a post-critical region where Va > i^arc- 

The real effect on the nucleation rate occurs when there is a transition 
from the pre-critical part of the source valley to the post-critical part of a 
destination valley. Transition from the post-critical part is useless because 
the embryos can simply continue to grow, they already overcame the barrier. 
So, there is no need to overcome another one barrier and this case is out of 
our interest. 
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At first we suppose that in the whole pre-critical part of the destination 
channel there is an equilibrium distribution. It means that there is no further 
change of channels and the destination channel will be the final destination 
channel. So, there is only one cascade - only one change of the channel. We 
shall call such processes as one-cascade processes. 

One can choose components in such a way that the first component is a 
rapid one. 

Consider the regions where the probability to change the channel is greater 
than to increase the value of slow components in the old channel. This cor- 
responds to condition 

Wsi<WiZ^e^v{Fr-F,) 

Here the kinetic coefficient Wsi is the total kinetic coefficient of all slow 
components, Wi is the kinetic coefficient of a rapid component and Zi is 
the corresponding Zel'dovich factor for transition over the ridge. The last 
inequality can be expressed in terms of the — F^ as 

Fr-F,< ln{Wsi/{WiZ^)) = At 

The rhs is a very slowly varying function. Approximately it is a constant. 
Consider 

A = Fr-F^ 

According to the approximate formulas the function A(z/a) has the second 
derivative 

which has a constant sign. 

Thus, A has no more than maximum (it will be marked by the index 
"m"). 

Certainly, the condition Fj. — Fy = At depends on the scale of Ua- It is 
necessary to choose the scale of fa-axis to have 

dA/dva ~ 1 

at A ~ At. Since A is not a too sharp function of Va-, it is easy to do. The 
condition A < A^ can be valid no more than in two zones: one before Uam 
another later Vam- Namely, in intervals 

[0,i/„t_], [vat+-,+oo\ 
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the last condition is valid. 

One can make the following notes: 

• The second interval [i'at+ , +00] can be absent when 

• The first interval also can be effectively (not precise) absent when 
Vat- < 1 which occurs rather often. 

• One can come to the situation when valleys are purely isolated. 

• One has to keep in mind that the approximate formulas take place only 
at big Ua. 

The interval [0, i/at_] will be called as the " pre-transition zone", the in- 
terval [i^at+y +00] - as the "post-transition" zone. 

Consider the question about the mutual position of the destination and 
the source channels. The definition of k as even without microscopic correc- 
tions {Y. ^4^1)7^'^^ contains 7 and Vi and is a very complex function. But in 
the majority of situations the increase of Va (other i/j are fixed) causes the 
increase of k. We shall imply this property to take place. This property will 
be referred as the property of K-convexity. The line k — const as a function 
of is convex. 

Certainly, in real systems there can be concave regions, where the growth 
of fa leads to the change of concentration, the partial volumes change, the 
surface tension change and the value of k falls. But this situation is exclusive. 

Under the property of convexity one can see that the transition will be 
carried out only by addition of molecules of the first component (ejection is 
not possible) and will go from the left side to the right side in Uf, plane. 

The precise position of boundaries have to be defined with surface ex- 
cesses. Also a shift connected with a special renormalization has to be taken 
into account. 

6.3 Nucleation conditions and supplying conditions 

Here we shall mark by x. 

Conditions for the possibility of nucleation through the post-transition 
zone are the following ones 



74 



• Transition has to be effective, i.e. there has to be a region in the 
post-critical region in a destination valley x > x^c where Ffi{x) < Fsc 
Certainly this region looks like 

[xfe, +00] 

and the beginning of this region has to be smaller than Xsc- 

Xf) <i ^sc 

• Transition has to be opened, i.e. 

The beginning of transition will be at 

There are two possibilities at Xy,: 

• The first possibility 

Here the common valley will be formed and the most effective transition 
will be at Xu defined as 

Fd{xu) = Fs{xy,) 

• The second possibihty 

Fdipi^w) ^ Fgixqj^) 

Here the transition from valley to valley occurs like a falling from the 
high channel to the low channel. Solution looks like the Trinkaus' one. 

To see the real process of the channel transition it is necessary to have 
corresponding conditions at the beginning of transition. These conditions 
have to be the equilibrium conditions. It is necessary that earlier in the 
valley there would be no possibility to escape from the valley. One has to 
analyze such possibility. 

To sec the transition in the pre-transition zone it is necessary that two 
conditions take place 
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• The first condition: 

Fd{x) < F,{x^) 

• The second condition: 

Transition has to lead to the post-critical region in the destination 
valley. 

We are interested to avoid such intensive transition which can destroy the 
equilibrium conditions at Xy^. 

Since has to be at the post-critical region, it is a decreasing function 
of X and it is sufficient to check condition at the boundary: 

Fd{xt-) < FsiXyj) 

In this situation the intensity of the valley transition in the pre-transition 
zone is so big that there is no equilibrium condition for the transition in the 
post-transition zone. 

Since the transition in the pre-transition zone has to lead to the post- 
critical zone then the peak of Fd lies inside the transition through the pre- 
transition zone. So, since 

maxFd > maxFs > Fs{xw) 

the transition occurs in a manner of common channel and the real transition 
takes place at Xp when 

Fd{xp) = Fs{xp) 

if 

Xp < Xf— 

If at 

X = Xt- 

we have 

Fd>Fs 

the most intensive transition takes place at Xt-. This situation is more prob- 
able than the precedent one. 

What has to be done when the condition 

Fd{xt-) > Fs{xt-) 
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takes place? 

Certainly, the transition can take place out of pre-transition and post- 
transition zones but with a very low probability. To take into account this 
possibility one has to add to Fg the quantity Fj. — Fg — Aj, i.e. to go from 
Fg to Fr — At. This has to be done out of pre- and post-transition zones. 

The point of transition will be near the root of equation 

Fa = Fr- At 

Let it be at x = Xy. 

This transition can not violate the equilibrium. So, the transition in the 
post-transition zone is not destroyed and intensities of this transition and 
transition in the post-transition zone have to be compared (added). 

We shall call this transition as "the saturation transition". 

Here the transition is going across the ridge into the valley. The surface 
excesses can be taken into account very simply by noticing that the forms of 
ridge and valley profiles remain the old ones and only the shifts of profiles as 
a whole take place due to the account of excesses. 



6.4 Details of the saturation transition 

Solution of the saturation transition is rather simple. Consider at first the 
general situation. Let nd{i'a) be the embryos number density in a destination 
valley, Ug^i'a) be the embryos number density in a source valley. The evolution 
equation for the source valley looks hke 

eM-Fr + Fs)W+ + exp(-F, + Fd)W+ 
Here Wa is kinetic coefficient, is the equihbrium distribution, the flow 

ng^eM-Fr + Fg)W+ 

is the flow from the source valley to the destination valley and 

n,|^exp(-F, + F,)Pr+ 
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is the flow from tlie destination valley to the source valley. The value Z is 
the Zeldovich' factor, A is the normalizing factor. 

Analogously one can write equation for the destination valley 



n,^ exp(-F, + - exp(-F, + F,) W+ 

In continuous approximation 



-n,^ exp(-F, + Fs)W+ + n,^ exp(-F, + Fd)W+ 



for the source valley and 



dud + dFd d 



rid-^ exp(-F, + + n,-^ exp(-F, + Fs)W^ 

^d 



for the destination valley. 
One can assume that 



dFd dFs 



(55) 



(56) 



are constants. Also it can be assumed that the linear approximations 

-Fr + Fd^ AdX + Cd (57) 

-Fr + Fs^ Asx + Cs (58) 

for X — Va ~ ^ao are valid. Here z/qq is some parameter chosen as to belong 
to effective region of transition. 

Then the stationary solutions will satisfy the system of equations 

^2 d 

ns + Vg-^ris - CgUs exp(^sx) + CdUd exp(^d2;) = (59) 



ox'^ ox 
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d 

g^^d + ^'^Q^^'^ ~ ^'^^'^ exp{Adx) + CsUs exp{Asx) = (60) 



with 



Cs = exp{Cs)ZsW,t/A,Was 

In the second solution because the region in the destination valley is the 
super-critical one it is possible to neglect -§^nd. Then the second equation 
is the linear first order differential equation with a known solution. Then 
after the substitution the first equation becomes the closed the closed equa- 
tion. Since solution of fl60l) contains the integral then to get the differential 
equation it will be necessary to differentiate fl59l) one time and the resulting 
differential equation will have the order 3. It can not be solved at least in 
elementary functions. So, it is necessary to consider simplification based on 
classification of transitions. 

These are three types of transitions - the non-equilibrium falling transition 
(first type), the equilibrium common valley transition (second type), the 
equilibrium saturation transition (third type). For different types we shall 
use different approximations. 

For the first type it is possible to neglect 

CdndfiXY>{Adx) 

in the first equation. Then it becomes the closed equation 
5^ d 

Us - Vs^ris - CsUs exp{Asx) = (61) 



ox'^ ox 

Then there is no necessity in validity of the linearization fl58|) . 

Solution of the last equation is presented above via cyhndrical function. 

The second equation is not necessary, but to complete the picture one 
can write it in the form 

d 

'"d.-Q^iT'd + CsUs exp(A^x) = (62) 



with the solution 



nd = - CsUs exp{Asx') /vddx' 
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Certainly, the presentation of the solution via the cylindrical function is 
not convenient. It is more convenient to fulfill a block transformation and 
then to solve the system of several algebraic equations. We shall formulate 
this procedure. Really, one can go from x to kx to have 

Ask ^ Ina 

where the parameter a ~ 1.5. Then with an increase of x by 1 the intensity 
of transition increases 1.5 times. Then the equation (16T!) will be 

d 

k~'^——ns + Vsk'^—Us - CsUso!^ = (63) 
ox^ ox 

Now it is possible to consider the interval —2 < x < 2 and to come back 
to the initial discrete form of equation 

k-^[nsix+l)-2ns{x)+ns{x-l)]+Vsk-^[ns{x+l)-ns{x-l)]/2-Csns{x)a'' = 

(64) 

These coupled algebraic equations have to be written at x = — 2,— 1,0,1,2. 
At X < — 2 one has to put to the equilibrium value. So, there is a system 
of five coupled equations which can be easily solved. 

One can continue analysis. Every band has a separate physical meaning: 

• The band x = 2 is the starting band. 

• At X = —1 one can use the smallness of the fiow Csns{x)a^ and the 
smallness of the deviation of from the equilibrium value. 

• The point x = is the point where \d/dx[{ns — n^)/Ti^]| attains maxi- 
mum and, thus, 

dydx^[ins-nl)/nl] = 

• At X = 1 one can assume that Us is already small in comparison with 
K 

• The values at x = 2 are the final values. 

One can use these features at the characteristic zones to get analytic solutions 
and then the common solution will be their combination. 

These approximations allow to solve this equation analytically by com- 
bination of the corresponding analytical band solutions. But the resulting 



80 



formulas will be very long to use them for calculations. To get concrete re- 
sults it is more profitable to solve algebraic equations, the precision is rather 
high while the error is less than one tenth. 

Certainly, one can directly solve the initial form of evolution equation 
described earlier as the discrete model. 

The analysis of the first type is completed. 

For the transition of the second type it is possible to neglect 

92 d 

ris + Vs—Us 



ox'^ ox 
and 



ox'^ ox 
Then both equations will be reduced to 

- CsUs exp(A^a;) + C7dnrfexp(Arfx) = (65) 

with the evident equilibrium solution as it was described earlier. The point 
a; = is the saddle point, i.e. the point where approximately Fg = Fa- 
The third type of transition can be described in a following manner: 



Equation (1591) leads to the fact that Ug = n^. For one can take 
approximation 

nlix) = nl{Q)exp{Bgx) 

where 



x=0 



dx 

Equation (I6OII looks like 
d 

^d-g^nd - CdJid exp(Arfa;) + C^n^ eyiY>{.^sx) = (66) 

and can be easily solved since it is the first order linear equation. The 
integral can be taken and the result will be expressed via Whittaker 
or Kummer functions which can be reduced to the Hypergeometric 
function. 
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The point x — has to be chosen as arg(maa;-J^). 

Since the result can be expressed only in the form of special functions it 
is worth solving the discrete model. The method is quite the same and one 
can come to 

_ <9^ _ 5 

k~'^Tr^nd + k ^va^r^d - C^rid exp(^d/cx) + Csns{0) exp((Ss + As)kx) = 
ox"^ ox 

(67) 

The value of k has to be chosen to satisfy 

k mm{Ad, As + Bg) — a 
Then algebraic equations will be 

k~^[nd{x + 1) - 2nd{x) + nd{x - 1)] + k~'^Vd[nd{xi) - nd{x - l)]/2 

(68) 

-Cdndexp{Adkx) + (^^^^(O) exp{{Bs + As)kx) = 

and have to be written at x = —2, —1, 0, 1, 2 

Also it is necessary to mention the possibihty to solve the discrete model 
from the very beginning. The starting equation will be 

W^[nd{x - 1) - "^^^^^^^Udix) - Udix) + ;^|f^^<^(^ + 1)] 

(69) 

-Cdndex.p{Adx) + Csnlex.p{Asx) = 

These equations are coupled algebraic equations. The initial condition is 
Ud — at X ^ — oo. 

Our consideration has to be completed by equation on parameters of 
transition. 

The points of approximations a; = form the equations on parameters of 
the process. The possible presence of the special functions can be eliminated 
by rational approximations for special functions. Then the parameters of th 
process will be determined by the root of the algebraic equation. 

Now one can see the general picture of transition. One can note that 
the complexity of the phase transition already between two valleys is rather 
essential. There are at least several possibilities to observe this transition 

• Equilibrium transition in the pre-transition zone 
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• Equilibrium transition out of transition zone 

• Equilibrium transition in the post-transition zone 

• Falling transition in the post-transition zone 

So, the unique approach to get the rate of nucleation is impossible. 

One has to stress that already the equilibrium transitions can lead to 
the absence of equilibrium in valley with bigger x, and n. This effect has 
to be taken into account to diminish the intensity of transition in the post- 
transition zone. 

Here it becomes clear that the flow can split and merge. Beside these 
effects one can see the rapid change of the leading manner of the supercritical 
embryo formation. This is caused only by kinetic coefficients and, thus, it is 
reasonable to speak about "the kinetic rupture in the rate of nucleation". 

One has to stress that in the saturation transition there is no difference 
whether the transition occurs in the open or in the forbidden zone. Really, 
Fj. — Fg can be greater than A^: 

Fr-Fs> At 

and the transition will take place. The only conditions is 

FriUa) - At < 

where F^c is the value at the saddle point with a minimal height and 

for some > Vacd 

Rigorously speaking the same consideration can take place for transition 

of other types. 

Then one can come to the situation when both the equilibrium common 
channel transition and the falling transition can take place. When at some 

> ^acd 

Fd{l^a) = Fs{Ua) < Fmc 

we have to examine F^. 
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If 

Fr — ^t< Fmc 

we see that the intensity of the common channel transition is greater than 
the intensity of transition through the saddle point. 

Since Fd for Ua > Vacd is a decreasing function this intensity is also greater 
than the intensity of a saturation transition. Then we have to compare it 
with the intensity of a falling transition. 

If 

- At < Fs 

we have the common valley transition with intensity greater than the further 
falling transition. 
If 

then the further falling transition will occur with intensity greater than the 
intensity of the common valley transition. 

All this is done without account of the soft shift. To take this shift into 
account in a rough approximation it is necessary to add to the quantity 

6.5 Other peculiarities of transition 

The property of the K-convexity is very important in the context of the 
current analysis because it forbids the possibility to reach the pre-critical 
region after the transition through the post-transition zone. Otherwise there 
would be a source of embryos in some region of the destination valley. The 
property of the K-convexity forbid the localization of the flow. 

Such a localization of the transition flow can be seen in a multi-valley 
transition. Consider the situation when there is an intermediate valley (index 
i) and, thus, there are two ridges - one between the source valley and the 
intermediate valley (index rs), another between the intermediate valley and 
the destination valley (index rd). Suppose that for intermediate valley 

Frs -Fi<At 

Frd -Fi<At 
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Then one can speak about one effective ridge with a height 

Then the property of the ridge convexity disappears and the locahzation of 
the transition flow can be seen. One can speak about the 

• Injection at the finite zone into the vaUey. 

Earher we consider only two components in the mixture. So, the inverse 
transition has to be the backward one. But in the three component mixture 
one can imagine the curved transitions - at first transition the concentration 
of the first component increases, at the second transition the concentration 
of the second component increases. However, it is necessary to have at least 
two rapid components in the mixture. In some special cases it is possible 
to return to the same valley but in another place of this valley - may be it 
is possible to jump from the pre-critical zone to the post-critical one, may 
be it is possible to make one loop of a spiral. Here the picture will be very 
picturesque. However, it would be very nice to see the concrete examples of 
such nucleating systems. 

Here we do not consider the transitions from the post-critical zone of one 
valley to the post-critical zone of the other valley because this transition can 
not change the rate of nucleation. 

Despite the transition will have now a very complex form the elementary 
bricks remain the same: 

• the equilibrium common-valley transition 

• the equilibrium saturation transition 

• the non-equilibrium falling transition 

The possibility to reach the rather transparent classification is based 
on the following simple approximate structure of every channel/separation 
line/ valley/ ridge: 

• Every channel/separation line/vallcy/ridge has a pre-critical zone which 
is directly (without hills) connected with the origin, post-critical region 
with the irreversible growth (until infinity) and a small near-critical 
growth. 
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The last property takes place for every channel, separation line, valley and 
ridge. 

According to the last analysis the multi-cascade transitions arc not ef- 
fective. Really, the cascade can lead to the post-critical region or to the 
pre-critical region. When it leads to the post-critical region it will be the 
last cascade. If it leads to the pre-critical region there is a smooth increasing 
path along the valley and this path will have at least the same intensity. So, 
the transition across the ridge is not effective here. As the result we see that 
there is only one main cascade in the multi-cascade transition. 

Here we imply that one cascade can be the saturation transition, the 
falling transition or the common valley transition. Actually, the saturation 
transition is also the common valley transformation because here there exists 
a common valley. Then we shall speak here about the generalized common 
valley. Then there is the generalized common valley and the falling transition. 

Certainly, the multi-common valley can be such a cascade. In this cascade 
at some may be finite zone several valleys are treated as one common channel. 
It is also possible that the set of common channels with the given channel 
can change - at some zone there is one set, at another zone there is another 
set. But in this common valley under the K-convexity there will be only one 
leading pair of the neighbor channels. 

As the result we see that in the binary case there is only one leading 
cascade which is the falling transition or the generalized common valley. 

We have examined only the stationary solutions. The relaxation of the 
distribution n to the stationary solution can be easily studied since in all sit- 
uations the stationary solution Ust is known. Then one can linearize equation 
on n — Ust and get 

dn 

where L is a differential operator (or in finite differences) on i/j. Then one 
can get the relaxation time as the minimal eigenvalue of the linear operator 
L in the evolution equation by the iteration procedures 

Trace{L), Trace{L'^), Trace{L^), etc. 

of the standard numerical methods. 

In reality all operators in the rhs of kinetic equations are reduced to 
the square approximations. Then the eigenvalues and eigenfunctions are 
known. Eigenfunctions are the Hermite's polynomials or the Generalized 
error-functions . 
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Main results 



One can see that the problem to get even the stationary rate of nucleation is 
rather complex. Below we present the sequence of actions to fulfill this task: 

1. We determine all channels and find the channel with a minimal activa- 
tion barrier. Determine its height Fcm 

2. We determine the rate Wa/Wb- Choose components to have Wa < Wb 
If Wb/Wa < exp{l) there will be a Stauffer's solution with Fcm- If there 
is an opposite situation one has to continue consideration. 

Suppose that we have the binary case and the ^-convexity. The last 
property is rather ordinary but it simplifies the consideration. Then 
the procedure will be the following 

(a) Instead of channels determine the valleys. D We determine also 
all ridges. We enumerate valleys to have < for i < j. We 
enumerate ridges to have < for i < j. For every neighbor 
valleys we determine the source valley i and the destination valley 
i + 1. Below we shall consider the onc-cascadc transition. 

(b) We determine the possibility of the saturation transition: there is 

satisfying conditions: 

^acs 

^ ^acr 
FriVa) — Af < Fcm 

If these conditions are satisfied we determine the point of the 
saturated transition t/* by equation 

This gives 

(c) We determine the possibility of the common valley transition: 
there exists Va with properties 
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J^a > ^acd 

The last condition determines only one point v** of the common 
valley transition with a maximal intensity. This value will be the 
saddle point of the unified valley. Here we determine 

If the equilibrium valley transition takes place there is no need to 
consider the falling transition. If it does not exist then we consider 
the falling transition. 

(d) The falling transition takes place when there is Va satisfying con- 
ditions: 

^acs 
> ^acd 
Fdil^a) < Fsi^a) 

FM - FM < At 

Conditions 

determine the point of transition v*** and the free energy 

F — F (u***) 

(e) To see what transition is more profitable one has to compare F^g, 
F^, F** and F^^^ and to choose the minimal value 

Fch — TTT'in'(^F(,g , Fjf, -F^**-^'**^,) 

This will be the free energy corresponding to this valley as the 
source valley. Then one has to take the minimum over all valleys 
and to determine the free energy of nucleation Ff. Then the rate 
of nucleation is rather approximately given by 

J = exp{Ff)ZWa 

where the Zeldovich' factor Z contains the normalizing factor of 
the equilibrium distribution. 
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(f) If there are two approximately equal minimal values of free en- 
ergies between Fes-, -^*) -^** and -F**^, then one has to add the 
quantity \n{d{Fr — Fs)/dv) to the free energy of the falling transi- 
tion. Certainly, the last quantity shifts the point of transition but 
approximately one can take it at the old point. 

The analysis presented above gave the following new results 

• The free energy of the embryo is found including the surface excesses 
and the clear interpretation of the generalized chemical potential is 
given. The variables giving the simple form of the free energy is found 
and their connection with the initial natural variables is shown (section 
1) . The similarity of the form of the near-critical energy to the situation 
without surface excesses is shown (section 1). The correction order of 
the Renninger-Wilemski's effect is shown (section 1). 

• The hierarchy in the near-critical region is shown (section 2). 

• The impossibility of the essential difference between the Reiss' formula 
ad the Stauffer's formula in capillary approximation is shown (section 
3). 

• The possibility to change valleys during the evolution is shown. The 
discrete analog of Trinkaus' solution is presented and investigated (sec- 
tion 4). 

• The possibility to have one united valley instead of several initial ones 
is shown. It is shown that the height of the effective activation barrier 
is changed in comparison with the heights of barriers in the initial 
channels. Thus, the rate of nucleation will be radically changed (section 
5) 

• The possibility to change the valley and to reach the post-critical zone 
already from the pre-critical transition is shown (section 6). This form 
a special type of the equilibrium saturation transition. This transition 
also leads to a new special value of effective height of activation barrier. 

• All possible transition are classified. It is shown that the tree mentioned 
types cover the variety of possible transitions. 
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Here only the main new results are outlined. An application of the pre- 
sented theory to the concrete binary and multicomponent systems will form 
the subject of a separate publication. 
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